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ABSTRACT 


Analytical  studies  have  been  made  of  laminar  film  condensation  on  a  horizontal 
elliptical  cylinder  in  a  pure  saturated  vapor  under  conditions  of  free  and  forced  convection. 
Estimation  of  interfacial  shear  stress  was  made  in  two  ways:  the  first  involving  an  asymptotic 
value  of  the  shear  stress  under  conditions  of  infinite  condensation  rate  and  the  second  based 
on  simultaneously  solving  the  two- phase  vapor  boundary  layer  and  condensate  equations.  The 
latter  approach  enables  the  determination  of  the  vapor  boundary  layer  separation  point.  For 
the  assumption  of  asymptotic  shear  stress,  effects  of  surface  tension  and  pressure  gradient  in 
the  condensatt  film  have  also  been  included.  At  the  extremes  of  eccentricity,  corresponding 
to  a  circular  tube  and  a  vertical  flat  plate,  the  results  are  compared  with  theoretical  and 
experimental  work  of  others.  Improvement  in  the  condensation  heat  transfer  coefficient  was 
found  for  elliptical  tubes  under  both  free  and  forced  convection  conditions  when  compared 
to  circular  tubes  of  the  same  surface  area.  In  the  latter  case,  this  improvement  was  due  mainly 
to  the  reduced  drag  of  the  elliptical  tube  providing  a  higher  vapor  velocity  for  the  same 
pressure  drop  as  that  across  a  circular  tube. 
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I.  INTRODUCTION 


A.  MOTIVATION  FOR  INVESTIGATION  OF  ELLIPTICAL  TUBES 

For  many  years,  condenser  design  in  the  United  States  has  been  based  on  the  Heat 
Exchange  Institute  (H.E.I.)  or  Tubular  Exchange  Manufacturers  Association  (T.E.M.A.) 
standards  [1,2]  which  rely  on  tube-side  empirical  relations  and  averaged  data  to  predict 
heat  transfer  parameters.  Since  there  is  a  large  degree  of  uncertainty  in  the  accuracy 
of  the  prediction,  the  condensers  are  over  designed  to  ensure  sufficient  margin  of 
reliability.  The  result  is  an  excessively  large,  expensive  condenser  for  the  desired 
thermal  duty  rating  requirements.  The  majority  of  condensers  currently  used  in  the 
Navy  were  designed  under  these  standards. 

Within  the  last  two  decades,  computer  modeling  has  been  used  to  more  accurately 
predict  the  heat-transfer  coefficients  of  condensers.  There  are  two  basic  approaches  to 
computer  modelling,  an  integrated  and  a  differential  approach,  of  which  the  former  is 
most  commonly  used  in  condenser  design  (Walker  [3]).  Integrated  approaches  divide  the 
condenser  into  zones  with  associated  mean  heat  transfer  and  fluid  properties  which  are 
typically  determined  from  semi-empirical  relations.  These  zones  are  integrated  over  the 
whole  condenser  and  provide  overall  performance  predictions  with  respect  to  shell  side 
pressure  drop  and  duty  rating.  Though  easy  to  use,  integrated  methods  are  restrictive 
when  used  in  conditions  which  are  outside  of  the  fluid  parameter  and  geometry 
constraints  of  the  empirical  relations  employed  in  the  approach.  Differential 

approaches  solve  the  basic  equations  of  fluid  mechanics  and  heat  transfer  and  provide 
an  understanding  of  local  variations  within  the  condenser,  but  are  too  complex  for  use 
in  condenser  design.  The  improvements  in  accuracy  as  a  rsuit  of  computer  design 
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methods  are  a  function  of  the  empirical  or  theoretical  models  used  in  their  development. 
Marto  [4}  provides  a  comprehensive  overview  of  the  state  of  computer  design  methods. 

The  purpose  of  empirical  and  theoretical  research  into  condenser  performance  is 
to  more  accurately  predict  condenser  characteristics  such  that  the  designer  can  develop 
reliable,  smaller,  and  hence  less  costly  condensers  which  meet  the  same  thermal  duty 
requirements.  Analyses  of  single  tubes  allow  isolation  of  individual  factors  which 
affect  condensation  without  the  complexities  which  result  from  vapor  flow  over  tube 
bundles  or  condensate  inundation.  Since  the  pioneering  work  of  Nusselt  [5],  a 
significant  amount  of  theoretical  work  has  been  done  analyzing  laminar  film 
condensation  on  horizontal  circular  tubes  (see  Rose  [6]),  the  results  of  which  have 
greatly  contributed  to  the  ability  to  reasonably  predict  single  tube  heat  transfer 
performance. 

Over  the  last  two  decades,  major  effort  has  been  expended  to  study  effects  which 
enhance  laminar  film  condensation  heat  transfer  over  that  obtained  from  a  plain 
circular  tube.  The  majority  of  these  techniques  have  focussed  on  controlling  the 
thickness  of  the  condensate  film  as  this  is  the  major  resistance  to  heat  transfer.  The 
condensate  film  thickness  generally  increases  with  streamwise  distance  from  the  top  of 
the  tube  and  is  dependent  on  both  the  rate  of  condensation  and  the  interfacial  shear 
between  the  vapor  boundary -layer  and  condensate  film.  Additionally,  at  high  vapor 
velocities,  vapor  boundary- layer  separation  occurs  which  results  in  a  rapid  thickening 
of  the  condensate  film  downstream  of  the  separation  point.  Enhancement  techniques 
include,  but  are  not  limited  to,  extended  surfaces  (fins),  profiled  tube  surfaces  (roped 
or  corrugated  tubes),  and  non -circular  tube  geometries.  Marto  [7],  Bergles  et  al.  [8]  and 
Webb  [9]  provide  comprehensive  reviews  on  such  enhancement  techniques.  These 
techniques  enhance  heat  transfer  through  an  increase  in  surface  area  to  volume  ratio 
and/or  the  use  of  surface  tension  to  thin  the  condensate  film. 
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Non-circular  geometries  contribute  to  a  thinning  of  the  condensate  film  by  placing 
more  of  the  surface  in  line  with  the  direction  of  gravity;  thus,  making  the  average 
streamwise  gravity  component  in  the  momentum  balance  of  the  condensate  film  larger 
than  that  obtained  from  a  circular  geometry.  Dhir  and  Lienhard  [10]  applied  a  simple 
Nusselt  type  analysis  to  an  arbitrary  plane  or  axisymmetric  body  in  which  the 
streamwise  gravity  component  varied  with  streamwise  length.  The  results  of  this 
analysis  gave  expressions  for  condensate  film  thickness  and  Nusselt  number  which  were 


identical  to  those  derived  by  Nusselt  except  that  the  gravitational  constant,  g,  is 


replaced  by  an  effective  constant,  geff,  which  is  defined  by 

-  x(g  H)* 

*  fJg'Wdx 


(1-1) 


where  x  is  the  streamwise  distance  from  the  leading  edge  of  the  surface  and  R  is  the 
radius  of  curvature  of  the  axisymmetric  body.  These  expressions  can  only  be  applied 
to  systems  for  which  curvilinear  coordinates  are  applicable,  i.e.,  for  bodies  in  which  the 
radii  of  curvature  are  much  greater  than  the  film  thickness. 

Shklover  and  co-workers  [11,12]  analyzed  a  horizontal  cylinder  with  a  logarithmic 
spiral  type  of  surface  curvature  (see  Figure  I-l)such  that  more  of  the  body  surface  is 
aligned  in  the  direction  of  gravity  and  the  radius  of  curvature  is  continuously 
increasing  along  the  streamwise  direction.  They  manufactured  this  tube  by  "mechanical 
deformation  of  a  circular  tube."  The  effect  of  this  profile  was  to  increase  the 
streamwise  effective  gravity  component  over  that  found  for  a  circular  cylinders  and 
to  utilize  surface  tension  effects  on  the  pressure  gradient  to  increase  the  mean  film 
velocity,  u(x,y),  in  the  streamwise  direction.  Their  analysis  considered  viscous,  gravity 
and  surface  tension  forces  in  the  momentum  balance  of  the  condensate  film  and 
neglected  vapor  shear  at  the  condensate/ vapor  interface.  These  effects  resulted  in  a 
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Figure  1-1.  Logarithmic  Spiral  Tube  from  Reference  (12]. 


thinning  of  the  condensate  film  and  an  increase  of  20-30%  in  the  overall  heat-transfer 
coefficient  when  compared  to  a  circular  tube  of  the  same  surface  area.  These  results 
agreed  with  experimental  data.  In  analyzing  the  individual  effects  of  gravity  and 
surface  tension,  it  was  determined  that  surface  tension  was  most  significant  over  the 
top  portion  of  the  tube  where  dR/dx  is  largest  and  accounted  for  10-15%  of  this 
increase. 

An  elliptical  tube  with  major  axis  aligned  with  the  direction  of  gravity  has  a 
larger  effective  gravity  than  a  circular  cylinder,  as  well  as  an  increasing  radius  of 
curvature  over  the  top  half  of  the  tube.  In  addition,  they  are  more  practical  from  a 
manufacturing  standpoint  than  the  non-circular  tubes  analyzed  above.  Vapor 
boundary -layer  separation  should  also  occur  at  a  point  further  downstream  as  compared 
to  a  circular  tube  and  thus  the  rapid  thickening  of  the  condensate  film  would  be 
delayed.  Methods  of  manufacturing  elliptical  tubes,  and  associated  problems  will  be 
discussed  in  Chapters  IV  and  VI. 


Wallis  [13]  conducted  flow  visualizations  (Figure  I-2)of  water  around  circular  and 


elliptical  tube  bundles.  These  photographs  provide  a  qualitative  feel  for  the  flow 


Figure  1-2.  Flow  Visualization,  Water  over  Circular  and  Elliptical  Tube  Banks, 
from  Reference  [13]. 


characteristics  around  each  type  of  tube  bundle.  In  examining  the  photographs,  it  can 
be  concluded  that  the  separation  point  for  the  elliptical  tube  bank  shifts  downstream 
compared  to  the  circular  tube  bank  and  results  in  smaller  wake  regions  behind 
individual  tubes.  These  effects  should  contribute  to  better  heat  transfer  performance 
for  the  elliptical  tube  bank.  Experimental  studies  by  Joyner  and  Palmer  [14]  show  that 
single  phase  flow  resistance  and  pressure  drop  are  significantly  smaller  for  elliptical 
tube  bundles  as  compared  to  circular  tube  bundles.  In  power  condensers,  shell-side 
(vapor)  pressure  drop  has  an  effect  on  the  overall  efficiency  of  the  power  plant  and  on 
the  thermal  driving  force  for  condensation  (TMI-TW1„).  Minimizing  the  vapor  pressure 
drop  improves  power  plant  efficiency  and  maintains  a  relatively  constant  Tsat  (for  a 
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large  pressure  drop,  Tsat  is  reduced  as  vapor  flows  through  the  condenser  resulting  in 
less  thermal  driving  force). 

Studies  involving  single-phase  heat  exchangers  and  condenser-evaporators  indicate 
that  elliptical  tube  heat  transfer  performance  is  superior  to  a  circular  tube  of 
comparable  surface  area.  Ota  and  Nishiyama  [15)  conducted  an  experimental  investi¬ 
gation  of  single  phase  forced  convection  heat  transfer  of  air  over  an  elliptical  cylinder 
of  minor-to-major  axis  ratio  of  1:3  at  various  angles  of  attack.  They  concluded  that 
elliptical  cylinders  gave  improved  heat  transfer  performance  over  circular  cylinders  at 
all  angles  of  attack  (as  a  result  of  fluid  turbulence)  in  the  range  of  Reynolds  numbers 
studied  (8000  -  79000).  Moalem  and  Sideman  (16)  conducted  a  theoretical  analysis  of  a 
horizontal  elliptical  tube  in  a  condenser-evaporator  used  in  desalinization  plants.  Their 
study  showed  a  10-20%  enhancement  for  elliptical  tubes  as  compared  to  circular  tubes 
with  the  maximum  enhancement  achieved  at  a  minor-to-major  axis  ratio  (aspect  ratio) 
of  1:4  where  the  major  axis  is  aligned  with  the  vertical.  Huang  and  Mayinger  [17] 
conducted  an  experimental  free  convection  heat  transfer  study  around  elliptical  tubes 
and  found  optimum  improvement  in  heat  exchanger  performance  for  tubes  with  major 
axis  aligned  vertically  and  with  an  aspect  ratio  of  1:2.  Merker  and  BShr  [18]  used  an 
analogy  between  momentum  and  heat  transport  to  derive  a  semi-empirical  relation  from 
which  to  determine  mass  transfer  rate  from  the  fluid  boundary-layer  into  the  free 
stream.  The  results  of  their  study  showed  improved  heat  transfer  performance  for  heat 
exchangers  in  which  the  elliptical  tubes  were  spaced  wider  in  the  longitudinal  direction 
and  more  compact  in  the  transverse  direction.  In  all  these  studies  it  appears  that  an 
elliptical  geometry  improves  the  heat  transfer  performance  of  a  heat  exchanger. 
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B.  SIGNIFICANT  THEORETICAL  STUDIES  ON  HORIZONTAL  CIRCULAR 


TUBES 

The  purpose  of  this  survey  is  to  provide  a  historical  perspective  to  the 
understanding  of  laminar  film  condensation  on  a  horizontal  circular  cylinder.  Only 
those  works  which  have  made  significant  advances  in  this  understanding  are  considered. 
The  methods  of  some  of  these  analyses  will  later  be  applied  to  formulate  and  solve  the 
governing  equations  for  the  case  of  the  horizontal  elliptical  cylinder.  Figure  1-3 
provides  the  geometrical  layout  and  coordinate  systems  used  in  these  analyses. 


Figure  1-3.  Geometry  for  Film  Condensation  on  a 
Horizontal  Circular  Cylinder. 


Nusselt  [5]  used  a  simple  momentum  and  energy  balance  to  determine  the  heat 
transfer  properties  for  condensation  on  flat  plates  and  horizontal  circular  cylinders. 
In  simplifying  these  physical  laws  for  solution,  he  made  the  following  assumptions: 


(1)  The  tube  wall  temperature  is  constant. 

(2)  The  thermophysical  properties  of  the  fluid  are  constant. 

(3)  The  temperature  at  the  film /vapor  interface  is  Tsat. 

(4)  The  condensate  film  thickness  is  small  compared  to  the  radius  of  the  tube. 

(5)  Heat  transfer  in  the  condensate  film  is  one  dimensional  in  the  radial  direction 
which  implies  a  linear  temperature  distribution  when  the  film  is  very  thin. 
Convection  in  the  condensate  film  is  neglected. 

(6)  The  only  forces  acting  on  the  condensate  element  are  due  to  viscosity  and 
gravity.  Inertial  forces  are  neglected. 
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(7)  Flow  in  the  condensate  film  is  laminar. 

(8)  The  vapor  is  quiescent  with  no  interfacial  shear. 


The  governing  equations  for  this  system  become 


— -  +  P i  8  sin<t>  =  0 
dx 1 


=  —(7V 


(1-2) 


where  m  is  the  condensate  mass  flux  rate  across  the  interface.  The  mean  heat  transfer 


is  given  by 


'•n* 


0.943 


kf  g  hk  p )  sin4> 


n,  L  A T 


1/4 


(I -3a) 


for  a  flat  plate  at  angle  4>  to  the  horizontal  (not  valid  at  <j>  =  0)  and 


<tffm  =  0.728 


Pi' 


il,  D  AT 


1/4 


(I-3b) 


for  a  horizontal  circular  cylinder. 

Sparrow  and  Gregg  [19,20]  extended  the  Nusselt  [5]  analysis  by  applying  boundary 
layer  theory  to  the  condensate  film.  This  study  incorporated  convection  and  inertia  in 
the  energy  and  momentum  balance  while  continuing  to  neglect  interfacial  shear.  They 
determined  that  for  practical  engineering  fluids,  inertia  and  convection  had  negligible 
effect. 

Shekriladze  and  Gomelauri  [21]  included  interfacial  shear  by  considering 
momentum  transfer  across  the  interface  due  to  the  condensation  process.  They  used  an 
asymptotic  expression  for  the  interfacial  shear  based  on  an  infinite  condensation  rate 
given  by 

t4  =  m(U *  -  u6 )  (I-J> 

where  is  the  streamwise  velocity  of  the  vapor  at  the  outer  edge  of  the  boundary- 
layer.  It  was  also  assumed  that  U^»Uj  such  that  the  film  velocity  may  be  neglected  and 
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that  inertia  and  convection  effects  in  the  condensate  film  could  also  be  neglected.  As 
noted  by  Rose  [6],  the  validity  of  neglecting  u5  requires  that  A.jAT/t^hfg  be  small  (large 
Prandtl  number).  For  forced  convection  heat  transfer  on  a  horizontal  circular  tube, 
Shekriladze  and  Gomelauri  determined  the  following  expression: 

Nu  =  0.9\[^e  (,'5) 

where 


Re  = 


P,  U.  D 

*1/ 


(1-6) 


is  the  two-phase  Reynolds  number.  When  the  effect  of  gravity  cannot  be  neglected 
(mixed  convection)  the  mean  heat-transfer  coefficient  is  given  by 

Nu  =  +V1  +  l-69  F  (l'?) 

where 


F  _  ^  hk  *  D 
X,  AT  Ul 


(1-8) 


and  measures  the  relative  effects  of  vapor  velocity  and  gravity.  It  should  be  noted  that 
based  on  potential  flow  outside  the  vapor  boundary  layer,  the  interfacial  shear  will 
always  be  positive  and  therefore  this  method  does  not  predict  vapor  separation  as  would 
occur  in  reality.  Based  on  a  minimum  separation  angle  of  82°  for  single  phase  flow  over 
a  circular  cylinder  without  suction  and  no  heat  transfer  after  separation,  Shekriladze 
and  Gomelauri  conservatively  estimated  that  the  heat-transfer  coefficient  would  be 
reduced  by  about  35%.  Actual  separation  occurs  at  a  point  further  downstream  and 
results  in  a  heat- transfer  coefficient  between  the  result  for  separation  at  82°  and  no 
separation.  These  expressions  have  been  found  to  be  reasonably  accurate  for  cases  of 
high  condensation  rate  and  are  useful  for  their  simplicity. 

Fujii  et  al.  [22]  conducted  an  analysis  of  mixed  convection  condensation  on 
horizontal  cylinders  using  the  two-phase  boundary-layer  equations.  Their  formulati  jn 
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neglected  inertia,  convection  and  pressure  gradient  effects  in  the  condensate  film.  The 
governing  equations  for  the  model  are  given  by 


dU 

_  4. 

3V  =0 

dx 

dy 

u *JL 

II 

+ 

dUk 

u* — - 

dx 

dy 

*  dx 

Pv  dy2 

(I -9a) 


for  the  vapor  boundary- layer  and 


du  ^  dv 
dx  dy 


0  =  h, — r  +  P,  g  sm<f) 
dy2 


X;AT 

TT 


fs 


(I -9b) 


for  the  condensate  film.  The  compatibility  equations  at  the  interface  are  given  by 

U*  =0 


An  approximate  integral  solution  of  the  momentum  equation  for  flow  over  a  cylinder, 
with  suction  by  Truckenbrodt  [23],  was  modified  to  more  closely  agree  with  the 
numerical  solutions  of  Terril  [24],  This  method  enabled  the  determination  of  the 
interfacial  shear  and  subsequent  solution  of  the  governing  equations  for  the  condensate 
film.  This  technique  also  predicted  the  point  of  vapor  separation  and  subsequently,  a 
more  accurate  heat- transfer  coefficient.  For  locations  downstream  of  the  vapor 
boundary- layer  separation  point,  the  interfacial  shear  was  assumed  to  be  negligible  and 
simple  Nusselt  theory  used  to  determine  the  heat  transfer.  Local  and  mean  Nusselt 
numbers  were  numerically  determined  for  boundary  conditions  of  uniform  wall 
temperature  and  uniform  wall  heat  flux.  Their  results  were  in  fair  agreement  with 
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with  experimental  data.  In  addition  to  the  dimensionless  parameter,  F,  a  new 
dimensionless  parameter,  G,  was  defined  as 

G  =  AT  k>  Pl  ^  (I- 10) 

n,A*\pvriy 

which  is  a  measure  of  the  condensation  rate.  As  G-°°,  the  results  of  Fujii  et  al.  [22] 
agreed  with  the  previous  results  of  Shekriladze  and  Gomeiauri  [21].  Lee  and  Rose  [25] 
noted  an  error  in  the  formulation  of  Fujii  et  al.  [22]  in  that  the  interfacial  shear  terms 
in  the  condensate  film  momentum  equation  should  have  been  divided  by  2.  They  re¬ 
analyzed  the  system  of  equations  using  both  the  modified  and  unmodified  method  of 
Truckenbrodt  [23]  and  found  little  difference  whether  one  used  the  modified  or 
unmodified  method  of  Truckenbrodt. 

To  date,  Gaddis  [26]  has  conducted  the  most  comprehensive  study  of  condensation 
on  horizontal  circular  tubes  using  the  two- phase  boundary -layer  equations.  His  analysis 
neglected  surface  tension  in  the  momentum  equation  as  well  as  viscous  dissipation  and 
pressure  in  the  energy  equation  of  the  condensate  film.  He  determined  that  inertia 
effects  were  negligible  for  most  media  (with  the  exception  of  liquid  metals)  and  that 
convection  effects  were  somewhat  significant  for  viscous  liquids  with  Pr  »  1,  similar 
to  the  conclusions  of  the  simpler  analyses  of  Sparrow  and  Gregg  [19,  20].  He  identified 
three  regimes  of  behavior  for  flowing  vapors.  For  low  vapor  Reynolds  numbers,  Rev, 
gravity  effects  dominate  the  heat  transfer  and  Nusselt’s  analysis  adequately  represents 
the  heat  transfer  behavior.  For  high  Rev,  condensation  was  shear  controlled  and  the 
Nusselt  number  is  proportional  to  the  square  root  of  Rey  as  found  by  Shekriladze  and 
Gomeiauri  [21].  The  intermediate  region,  in  which  gravity  and  vapor  shear  are  both 
significant,  results  in  no  simple  relation.  Gaddis  also  analyzed  several  cases  of  flow 
separation  (vapor  boundary- layer  and/or  condensate  film  separation).  For  low  Rey  (  * 
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102),  vapor  boundary-layer  separation  occurred  with  no  flow  reversal  in  the  liquid  film. 
This  result  is  expected  as  the  pressure  gradient  effect  on  the  condensate  film,  which 
tends  to  slow  the  condensate  flow  over  the  rear  half  of  the  tube,  is  a  function  of  the 
vapor  potential  velocity.  At  moderate  Rev  (=  103)  vapor  boundary -layer  separation  and 
conder  .ate  flow  reversal  near  the  interface  occurred  with  a  resulting  decrease  in  the 
mean  fluid  velocity  and  rapid  thickening  of  the  condensate  film.  At  high  Rev  (=  105), 
no  vapor  boundary-layer  separation  occurred  but  the  condensate  film  separated  at  the 
wall  causing  a  sharp  reduction  in  mean  condensate  velocity  with  a  rapid  thickening  of 
the  film  as  a  result  of  a  significant,  adverse  pressure  gradient.  Fluid  separation  at  their 
respective  boundaries  is  defined  as  the  condition  at  which  the  shear  is  less  than  or  equal 
to  zero. 

diMarzo  and  Casarella  [27]  conducted  a  similar  analysis  using  a  more  general 
solution  technique  and  arrived  at  the  same  results  as  Gaddis  [26],  In  their  analysis  of 
the  flow  separation  phenomenon,  they  provided  guidance  on  determining  heat  transfer 
performance  once  separation  had  occurred.  At  low  Rev,  gravity  effects  dominate  and 
no  flow  reversals  or  thickening  of  the  condensate  film  occurred.  The  use  of  a  Nusselt 
type  analysis  with  no  interfacial  shear  would  seem  prudent  in  the  region  beyond  vapor 
separation.  At  moderate  to  large  Rev,  pressure  gradient  effects  dominate.  The  adverse 
pressure  gradient  over  the  back  half  of  the  tube  causes  flow  reversal  and  a  rapid 
thickening  of  the  condensate  film.  In  this  case,  it  would  be  prudent  to  neglect  heat 
transfer  completely  beyond  the  separation  point  of  the  vapor- boundary  layer  or 
condensate  film. 

Rose  [28]  studied  the  effects  of  pressure  gradient  in  the  condensate  film.  In  this 
case,  the  pressure  gradient  on  the  condensate  film  is  due  to  the  pressure  gradient  of  the 
vapor,  as  determined  by  potential  theory,  which  is  impressed  on  the  condensate  film. 
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He  simplified  the  analysis  by  using  the  Shekriladze- Gomelauri  model  (21]  which  ignores 
inertia  and  convection  in  the  condensate  film  and  uses  the  asymptotic  expression  for 
the  interfacial  shear.  As  a  result  of  this  analysis,  a  dimensionless  parameter,  P,  was 
defined  by 

P  =  (I- 11) 

p,  K 

which  represents  the  pressure  gradient  effect  (Gaddis  (26)  had  a  similar  combination 
of  dimensionless  parameters  which  were  equivalent  to  P).  When  P  =  0,  the  governing 
equations  reduce  to  those  of  Shekriladze  and  Gomelauri  (21).  He  concluded  that  the 
effect  of  including  pressure  gradient  was  to  improve  heat  transfer  over  the  forward 
half  of  the  tube  since  the  pressure  gradient  over  this  region  is  favorable  (tends  to 
increase  the  mean  film  velocity  which  reduces  the  condensate  film  thickness).  As  a 
result  of  the  formulation,  he  found  that  for  cases  where  P  >  F /8,  a  critical  angle  at 
some  point  on  the  rear  of  the  tube  was  reached  where  d&‘/d<t>  -  <*>.  In  this  case,  it  was 
not  possible  to  obtain  a  solution  over  the  entire  tube.  It  was  postulated  that  this  critical 
angle  might  indicate  some  instability  followed  by  some  degree  of  waviness  or 
turbulence.  It  was  noted  that  this  critical  point  was  reached  prior  to  the  point  at  which 
the  condensate  film  separated.  For  conditions  which  permitted  solution  over  the  entire 
tube  (i.e.  P  <  FI 8),  any  increase  in  heat  transfer  achieved  over  the  forward  part  of  the 
tube  was  balanced  by  a  decrease  in  heat  transfer  as  a  result  of  the  adverse  pressure 
gradient  over  the  back  half  of  the  tube  such  that  there  was  little  change  in  the  mean 
heat-  transfer  coefficient.  Since  the  pressure  gradient  due  to  potential  flow  of  vapor 
just  outside  the  vapor  boundary- layer  is  given  by 


dP 

dx 


-p/V 


dH± 

dx 


(M2) 
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which  is  symmetric  around  the  surface  of  the  tube  from  front  to  back,  the  favorable 
effect  over  the  front  half  of  the  tube  is  exactly  offset  by  the  adverse  effect  over  the 
back  half  of  the  tube.  Referring  to  the  studies  conducted  by  Gaddis  [26J  and  diMarzo- 
Casarella  [27],  in  the  region  of  moderate  Rev,  condensate  film  separation  from  the  tube 
wall  does  not  occur  but  the  pressure  gradient  sufficiently  retards  the  film  such  that 
flow  reversals  occur  and  the  film  still  rapidly  thickens.  In  the  present  study,  the  effects 
of  pressure  gradient  using  the  analysis  of  Rose  [28]  resulted  in  a  condensate  velocity 
distribution  which  approached  zero  as  the  critical  angle  was  reached.  This  phenomenon 
may  be  a  result  of  the  model  used  in  formulating  the  problem.  Since  Shekriladze  and 
Gomelauri  [21]  used  an  asymptotic  expression  for  the  interfaciai  shear,  negative 
velocities  cannot  exist  in  the  velocity  profile.  The  positive  velocity  distribution  is 
retarded  by  the  adverse  pressure  gradient.  In  the  case  of  P  >  F 1 8,  the  pressure  gradient 
is  large  enough  to  cause  complete  stoppage  of  film  flow  at  <J»C  and  the  film  thickness 
increases  rapidly. 

Krupiczka  [29]  examined  the  effects  of  surface  tension  due  to  film  curvature  on 
condensation  on  circular  cylinders.  He  used  a  simple  Nusselt  type  model  but  included 
surface  tension  in  the  momentum  equation.  In  his  development  he  did  not  assume  the 
film  thickness  was  much  less  than  the  radius  of  the  cylinder  until  after  the  inclusion 
of  the  surface  tension  term  to  account  for  the  curvature  of  the  film.  The  resulting 
equation  was  a  second  order  ordinary  differential  equation  which  required  two  initial 
conditions.  The  first  condition  was  given  by  the  symmetry  of  the  problem.  However, 
the  initial  thickness  was  not  obvious  and  was  arbitrarily  chosen  to  be  that  obtained 
from  the  Nusselt  model.  He  concluded  that  the  effect  of  surface  tension  was  small  on 
the  forward  part  of  the  tube  and  increases  in  significance  over  the  k  part  of  the  tube 
due  to  the  rapidly  changing  film  thickness.  This  significance  was  dependent  on  the 
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magnitude  of  dimensionless  parameter,  A,  given  by 

A  =  -We 
4 

which  is  dependent  on  the  surface  tension  of  the  fluid  through  a  "modified”  Weber 
number  (We  -  o/( pgr2))  and  the  radius  of  the  cylinder,  r.  The  conclusion  was  that 
surface  tension  can  be  significant  for  small  diameter  tubes  and  wires  where,  as  r 
becomes  smaller,  the  film  thickness  becomes  relatively  significant  and  hence  the  film 
curvature  cannot  be  neglected.  For  most  practical  cases,  and  therefore  the  surface 
tension  effect  due  to  film  curvature  can  be  neglected.  The  effect  of  assuming  the 
initial  film  thickness  to  be  that  given  by  Nusselt  [5]  was  checked  and  verified  to  have 
minimal  effect  on  the  mean  heat- transfer  coefficient. 

C.  SIGNIFICANT  THEORETICAL  STUDIES  ON  HORIZONTAL  ELLIPTICAL 

TUBES 

Some  theoretical  work  has  been  done  on  horizontal  elliptical  cylinders  in  a 
quiescent  pure  vapor  using  a  Nusselt  type  model.  The  physical  orientation  and 
coordinate  system  is  shown  in  Figure  1-4. 

Cheng  and  Tao  [30]  approximated  the  surface  of  an  ellipse  by  several  circular  arcs. 
They  analyzed  condensation  on  these  arcs  using  the  same  assumptions  as  Nusselt.  The 
ellipse  was  aligned  such  that  the  major  axis,  a,  was  aligned  with  the  direction  of 
gravity.  From  their  numerical  results  they  determined  that  the  heat-transfer 
coefficient  decreased  with  increasing  eccentricity,  k  (defined  as  the  ratio  of  the  minor 
to  major  axis).  Values  of  k  were  varied  from  that  of  a  vertical  flat  plate  (k  =  0)  to  that 
for  a  horizontal  circular  cylinder  (k  =  1).  The  mean  Nusselt  number  was  determined  by 
a  surface  area  weighted  average  of  the  mean  Nusselt  number  for  each  circular  arc  and 
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Figure  1-4.  Geometry  for  Film  Condensation  on  a  Horizontal  Elliptical  Cylinder, 
is  given  by 

Nu  =  0.2876  [—  V4  £  |-P  ■  ~~  g-  f9uZ,(6 VAdQ  (I  *  I4) 

Ja\  [  ti? 

where 

Zfs J»)  =  -(sine^ffsmO)^  ♦  Cf  .  (I- 15) 

3  J8U 

Cj  is  an  integral  constant  which  links  each  arc  and  i  is  the  number  of  arcs  used  in  the 
approximation.  In  the  practical  range  of  eccentricity  (0.3  -  0.6),  the  mean  heat  transfer 
coefficient  was  increased  by  10  to  18%  over  that  of  a  circular  tube  with  the  same 
surface  area.  Ali  and  McDonald  [31]  conducted  a  similar  type  analysis  as  Cheng  and 
Tao  without  the  circular  arc  approximation  as  a  first  estimate  for  condensation  on 
inclined  circular  tubes. 

Wang  et  al.  [32]  used  the  Nusselt  assumptions  to  analyze  condensation  on 
horizontal  elliptical  tube  for  which  the  major  axis  is  oriented  at  an  angle,  a.  with 
respect  to  the  vertical  axis.  They  obtained  an  expression  for  the  mean  heat- transfer 
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coefficient  given  by 


where 


a 


1  6ft  P  ]  8 
4  x],  AT 


1 

yjab 


u* 

d*  *  r«»>w  d») 


(1-16) 


Zi  P)  =  sin^p  f*  I 


T3/2 


a2sin2(P-a)  +  f>2sin2(P~a) 


sini/3p  dp 


(1-17) 


The  results  of  this  study  showed  that  the  maximum  mean  heat- transfer  coefficient  is 
obtained  when  the  major  axis  is  aligned  with  the  vertical  (a  =  0°  ).  These  theoretical 
results  were  validated  with  experimental  data  for  an  elliptical  tube  with  semi -minor 
axis  of  8  mm  and  semi- major  axis  of  22  mm.  Sheng  and  Cha’o  [33]  noted  that  the  mean 
heat- transfer  coefficient  was  incorrectly  determined  since  the  heat- transfer  coefficient 
is  averaged  over  the  surface  area  and  the  radial  distance,  r,  is  not  constant  for  an 
ellipse.  In  the  course  of  this  current  study,  the  problem  formulation  of  Wang  et  al.  [32] 
was  rerun  with  the  correct  expression  for  the  mean  heat-transfer  coefficient  (as  later 
derived  in  Chapter  II).  The  corrected  mean  heat- transfer  coefficient  of  the  elliptical 
tube  analyzed  by  Wang  et  al.  is  11.3%  larger  then  a  circular  tube  with  the  same  surface 
area.  Additionally,  it  should  be  noted  that  the  heat  transfer  performance  of  this 
elliptical  cylinder  was  better  than  the  circular  cylinder  for  angular  orientations  up  to 
a  =  50°.  The  corrected  theoretical  results  are  approximately  5%  lower  than  those 
determined  by  Wang  et  al.  and  more  closely  agree  with  their  experimental  data.  Figure 
(1-5)  shows  the  heat  transfer  enhancement  for  elliptical  tubes  of  varying  eccentricity 
and  orientation  angle. 

Sheng  and  Cha’o  [33]  studied  the  effects  of  surface  tension  and  variable  wall 
temperature  on  condensation  on  a  horizontal  elliptical  tube.  The  remaining  assumptions 
of  the  analysis  were  the  same  as  those  above.  For  the  wall  temperature,  a  cosine 
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Figure  (1-3).  Elliptical  Tube  Mean  Nu  Compared  to  Circular  Tube  Mean  Nu  for 
varying  Eccentricity  and  Orientation  Angle. 


distribution  of  the  form 

AT  =  -  A  cos0)  (M8) 

was  used  (which  has  been  shown  (Memory  and  Rose  [34])  to  be  in  good  agreement  with 
experimental  data)  where  the  value  of  A  depends  on  the  ratio  of  the  outside  to  inside 
heat-transfer  coefficient.  They  determined  that  variable  wall  temperature  affected 
local  values  but  not  the  mean  values  of  the  heat- transfer  coefficient.  For  the  surface 
tension  effect,  it  was  assumed  that  the  film  thickness  was  much  smaller  than  the  radius 
of  curvature  of  the  elliptical  surface  and  thus  the  surface  tension  effect  was  due  solely 
to  the  curvature  of  the  tube  wall.  This  is  in  contrast  to  the  study  of  Krupiczka  [29]  for 
a  circular  tube  where  the  surface  tension  was  due  to  the  curvature  of  the  condensate 
film.  Over  the  back  half  of  the  tube  the  radius  of  curvature  decreases  with  streamwise 
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distance  resulting  in  a  retarding  effect  on  film  flow.  As  the  flow  is  slowed,  the  film 
thickness  increases  rapidly  and  at  some  critical  angle  becomes  infinite  in  magnitude. 
This  phenomenon  is  similar  to  that  described  by  Rose  [28]  for  the  pressure  gradient 
effect.  The  surface  tension  effect  was  determined  to  be  significant  for  k  <  0.6  .  In 
comparison  to  a  circular  tube,  surface  tension  causes  a  favorable  pressure  gradient  over 
the  front  half  of  the  tube  resulting  in  a  greater  film  velocity,  thinner  film  thickness  and 
improved  heat  transfer.  It  has  an  opposite  effect  on  the  back  half  of  the  tube  which 
tends  to  negates  this  improvement.  The  gravity  component  of  the  momentum  balance 
is  the  driving  force  behind  the  enhancement  in  elliptical  tubes.  Consideration  of 
surface  tension  results  in  a  slight  decrease  in  the  mean  Nusselt  number  as  compared  to 
the  situation  in  which  the  surface  tension  is  neglected. 

A  potential  advantage  of  an  elliptical  tube  compared  to  a  circular  tube  is  the 
difference  in  vapor  flow  characteristics  as  a  result  of  a  better  streamlined  shape. 
Panday  [35]  developed  an  explicit  numerical  method  for  two  dimensional  film 
condensation  and  applied  it  to  the  case  of  downward  flowing  vapor  over  elliptical 
cylinders.  Convection  and  inertia  were  included  in  the  condensate  film  as  well  as 
surface  tension  and  pressure  gradient  (as  a  result  of  potential  flow  of  vapor  outside  the 
vapor  boundary- layer).  The  interfacial  shear  was  approximated  using  the  asymptotic 
expression  for  infinite  condensation  rate.  Several  errors  were  found  in  the  expressions 
for  surface  tension  and  pressure  gradient  as  a  result  of  an  incorrect  analysis  of  the 
differential  stream  wise  length,  dx.  The  first  error  is  the  result  of  assuming  that  dx  = 
rdtf>-  This  relationship  assumes  that  the  radial  distance  from  the  centroid  of  the  ellipse 
is  constant  over  the  interval  of  the  parametric  angle.  This  fact  is  not  true  (as  will  be 
shown  in  Chapter  II)  and  results  in  an  error  whose  magnitude  is  dependent  on  the  step 
size  used  in  the  numerical  procedure.  The  second  error  involves  the  expression  for  r. 
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If  the  parametric  angle  is  measured  from  the  vertical  axis  of  the  ellipse  (as  appears  to 
be  the  case  in  Panday’s  expression  for  vapor  velocity  potential),  then  r  should  be  given 
by 

x  =  b  sin$ 

y  =  a  cos<|>  (I  - 18) 

r  =  'Jx1  +  y2  =  n^cos24>  +  k2  sin2$ 

Panday  used  an  expression  for  radial  distance  given  by 

r  =  aysin2^  +  Jt2cos24>  (l- 19) 

The  magnitude  of  the  resulting  error  is  dependent  on  the  eccentricity  of  the  ellipse. 
Additionally,  there  appears  to  be  an  error  in  the  gravity  component  of  the  momentum 
equation.  Panday  uses  a  body  force  given  by  pgsinB,  which  is  true  only  for  circular 
geometries  but  is  untrue  for  elliptical  geometries.  His  results  will  therefore  not  account 
for  the  improved  performance  which  are  obtained  by  placing  more  of  the  surface  in 
line  with  the  direction  of  gravity.  Panday’s  conclusions  were  that  the  overall  heat 
transfer  was  reduced  for  elliptical  tubes  at  low  velocities  (as  compared  to  a  similar 
circular  tube)  due  to  a  rapid  thickening  of  the  film  at  the  stagnation  point  and  the 
overall  heat  transfer  is  increased  for  high  velocities  due  to  increased  interfacial  shear. 
These  conclusions  are  opposite  to  what  one  may  expect.  At  low  velocities,  the  film 
thickness  near  the  stagnation  point  for  elliptical  cylinders  is  thinner  than  the  equivalent 
circular  tube  as  a  result  of  increased  effective  gravity.  This  effect  results  in  improved 
heat  transfer.  Pressure  gradient  effects  due  to  potential  flow  and  interfacial  shear 
should  be  negligible  since  the  velocity  is  small.  At  high  velocities,  the  streamlined 
geometry  of  elliptical  cylinders  results  in  lower  interfacial  shear  and  pressure  gradient 
effects  as  compared  to  circular  cylinders.  These  effects  should  result  in  a  thickening 
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of  the  condensate  film  as  compared  to  circular  cylinders  and  hence  reduced  heat 
transfer  for  the  same  effective  diameter  and  free  stream  vapor  velocity. 

D.  OBJECTIVES  OF  THIS  THEORETICAL  STUDY 

The  results  of  previous  studies  indicate  that  elliptical  tubes  can  be  used  to 
increase  condensation  heat  transfer  in  low  vapor  velocity  condensers  as  compared  to 
circular  tubes  with  the  same  surface  area.  The  present  investigation  examines  the 
effects  of  vapor  shear,  pressure  gradient  and  surface  tension  in  laminar  film 
condensation  on  a  single  horizontal  elliptical  tube  with  its  major  axis  aligned  with 
gravity  and  the  free  stream  velocity.  Interfacial  shear  is  estimated  using  both  the 
simple  assumptions  of  Shekriiadze  and  Gomelauri  [21]  and  the  more  complex  technique 
of  Fujii  et  al.  [22],  The  latter  case  calculates  the  angle  at  which  the  vapor  boundary- 
layer  separates  such  that  the  effect  of  reduced  drag  on  the  mean  heat-transfer 
coefficient  can  be  evaluated. 
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II.  THEORETICAL  DEVELOPMENT 

Following  the  theoretical  developments  used  for  the  case  of  laminar  film 
condensation  on  a  horizontal  circular  tube,  theoretical  models  are  developed  for  laminar 
film  condensation  on  a  horizontal  elliptical  tube  of  eccentricity,  k.  This  development 
starts  with  a  Nusselt  [5]  type  model  in  a  quiescent  vapor,  and  then  adds  forced  and 
mixed  convection  as  in  the  models  of  Shedkriladze  and  Gomelauri  [21]  and  pressure 
gradient  as  considered  by  Rose  [28].  The  pre-^ure  gradient  takes  into  account  the 
effects  of  potential  flow  outside  the  vap.  -oundary-layer  as  well  as  surface  tension. 
Finally,  a  model  is  developed  which  analyzes  the  vapor  boundary- layer  and  boundary- 
layer  separation  following  Fujii  et  al.  [22].  Pressure  gradient  and  surface  tension  are 
not  considered  in  this  model  due  to  the  complexities  introduced  by  the  use  of  the  two- 
phase  boundary- layer  equation.  Where  possible  these  elliptical  models  are  checked 
against  existing  theories  for  the  "limiting"  eccentricities  of  a  circular  cylinder  ( k=0 ), 
vertical  flat  plate  (k=l)  and  horizontal  flat  plate  (k  -<*>). 

A.  FACTORS  RELATED  TO  ELLIPSE  GEOMETRY 

Consider  an  elliptical  cylinder  whose  cross-section  is  oriented  such  that  the  major- 
axis  is  aligned  with  the  vertical  as  shown  in  Figure  II- 1.  The  eccentricity,  k,  of  the 
ellipse  is  defined  by 

*  =  £  .  (II-l) 

a 

where  a  and  b  are  the  semi-major  and  semi- minor  lengths,  respectively.  Functions 
related  to  the  geometry  of  the  ellipse  are  initially  developed  in  a  cartesian  coordinate 
system  (x^yj)  whose  origin  coincides  with  the  centroid  of  the  ellipse.  This  is  then 
transformed  into  a  cylindrical  coordinate  system  (r,6)  where  r  is  the  radial  distance 
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from  the  ellipse  centroid  to  a  point  on  the  ellipse  surface  and  0  is  the  angle  measured 


from  the  vertical.  The  transformation  equations  are 

x,  =  r  sin 0 

y ,  =  r  cos0  . 


(«- 2) 


The  elliptical  surface  can  also  be  defined  by  a  parametric  angle,  <j>,  measured  from 
the  upper  semi-major  axis,  such  that  the  Vj  coordinate  on  a  circle  of  radius  a  translates 


Figure  II- 1.  Geometry  for  Film  Condensation  on  a  Horizontal  Elliptical 
Cylinder  with  Major  Axis  Parallel  to  Gravity. 

to  the  yt  coordinate  on  the  ellipse  and  the  coordinate  on  a  circle  of  radius  b  translates 

to  the  Xj  coordinate  on  the  ellipse.  The  transformation  equations  are  given  by 

x,  =  b  sin<|> 

I  T  /IT  'l'V 


(II- 3) 


y,  =  a  cos4>  • 

Combining  Equations  (1-2)  and  (1-3)  results  in  a  relationship  between  r.0  and  <f>: 
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sind>  =  —  sin0 
b 

COS<t>  =  —  COS0 
a 


tan<t>  = 


tan0 


(a) 

(b) 

(c) 


(II-4) 


It  is  assumed  that  the  condensate  thickness,  6,  is  much  smaller  than  the  radius  of 
the  elliptical  surface.  Therefore,  the  equations  for  the  condensate  film  and  vapor 
boundary -layer  are  developed  using  an  orthogonal,  curvilinear  coordinate  system,  (x,y), 
where  x  is  the  streamwise  distance  along  the  elliptical  surface  and  y  is  the  distance 
normal  to  the  surface. 

1.  Radial  Distance,  r 


In  cartesian  coordinates,  the  surface  of  an  ellipse  is  given  by 

b2  a2 


(H-5) 


The  radial  distance,  r,  can  be  determined  by 


**!»?, )  =  /  A  +  y\ 


(H-6) 


Using  the  n-rametric  transformation  equations  (Equation  (II - 3)),  an  equivalent 
expression  for  the  radial  distance  is  given  in  terms  of  4>  by: 

K40  =  a\j  cos2^  +  k2  sin24>_  .  C11'7) 
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2. 


Radius  of  Curvature  of  an  Ellipse,  R 

The  curvature,  k(Xj),  and  the  radius  of  curvature,  R(X]),  of  a  function,  f(Xj), 


is  given  by 


*(*,)  = 


1 A*,)  I 


[1  +  (f'(x1))2]3'2 
R(x{)  =  1/k(Xj) 


(II  *8) 


where  the  prime  denotes  differentiation  with  respect  to  xt.  Solving  for  as  a  function 
of  Xj  using  Equation  (II -  5),  determining  the  radius  of  curvature  using  Equation  (II - 8) 
and  transforming  the  results  to  parametric  coordinates  gives  the  radius  of  curvature  as 
a  function  of  <|>: 

R(<t>)  =  ^  [sin2<|>  +  >t 2cos2<t>]3/2  .  (II -9) 


3.  Stream  wise  Length,  x 

Since  the  radial  distance  to  a  point  on  the  elliptical  surface  is  not  constant, 
the  streamwise  length  is  not  proportional  to  0  as  it  would  be  for  a  circle.  Additionally, 
as  eccentricity  varies,  the  streamwise  length  for  a  given  0  is  not  the  same.  To  obtain  an 
expression  for  x  in  terms  of  parametric  angle  <J>,  consider  first  a  point  on  the  ellipse 
surface  as  shown  in  Figure  II- 2.  Moving  a  small  distance,  dx,  results  in  incremental 

changes,  d0and  dr.  The  resulting  relationship  between  x,r  and  0  is  given  by. 

(dx)2  =  (dr)2  +  (r  d0)2  (II- 10) 

Expressions  are  now  needed  for  dr  and  dd  in  terms  of  d<f».  For  dr,  this  is  simply 

achieved  by  differentiating  Equation  (II - 7): 

dr  =  iJLgLjg  M”2*  dd>  .  (II-H) 

^  y/cos2^  +  Jt2sin24> 
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Figure  II-2.  Enlarged  View  of  Point  on 
Ellipse. 


For  d0,  taking  Equation  (ll-4a),  substituting  Equation  (II - 7)  for  r  and  differentiating 


tbe  resultant  expression  results  in: 


art* - - 

[cos2^  +  fc2sm24>]3/i 


(II- 12) 


Combining  Equations  (II -4b)  and  (II- 7),  and  substituting  for  cos0  in  Equation  (II- 12) 


results  in  an  expression  for  dd  in  terms  of  d<fr. 


dB  =  - - - d< {>  . 

cos2^  +  it2siii24> 


(II- 13) 


Substituting  Equations  (II- 11)  and  (11-13)  into  Equation  (II- 10)  gives  an  expression  for 


dx  as  a  function  of  d<fy. 


*2+-(ifc2-l)2siir2$ 

dx  -  a  - * - d$ 

\  cos2^  +  /k2sin24> 


(II- 14) 


This  may  be  integrated  to  obtain  the  streamwise  length.  For  compactness,  a  function 


X(4>)  is  defined  by 


X(4>)  =  a 


k2  +  — ()t2-l)2sm22$ 
_ 4 _ 

cos2^  +  k2sirr<|> 


(II- 15) 
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In  later  developments,  a  characteristic  length  is  used  to  non-dimensionalize 
the  heat  transfer  parameters  in  the  model.  An  effective  diameter,  De,  is  defined  as  the 
diameter  of  a  circular  cylinder  having  the  same  surface  area  as  an  elliptical  cylinder 


and  will  be  used  as  a  characteristic  length.  The  effective  diameter  is  given  by 


nDt  =  2  /fl"x(4>>*4>  ■ 


(11-16) 


Streamwise  length  is  non-dimensionalized  by: 


(11-17) 


This  enables  direct  comparison  between  elliptical  and  circular  cylinders  since  it 
represents  an  equivalent  surface  area. 

4.  Streamwise  Gravity  Component,  gz 

The  component  of  gravity  in  the  streamwise  direction  is  tangent  to  the 
elliptical  surface  and  is  a  function  of  streamwise  location.  A  line  tangent  to  the 


surface,  defined  by  the  slope  of  the  surface,  is  obtained  from  Equation  (II -5), 

m  a  *i 

A  unit  vector  in  the  streamwise  direction  is  then  given  by 


(11-18) 


(b2  -  x,2)  +  ^-x\ 


f7i  2  ^  a  ? 

V*  "  1  ‘  ^  *i  J 


(IM9) 


The  streamwise  gravity  component  is  given  by  the  dot  product  of  the  gravity  vector  and 


the  tangent  vector, 


8,  s  | 


(»!  -  -,1)  *  v? 

b 2 


(II -20) 
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Transforming  Equation  (II-20)  into  parametric  coordinates  results  in  the  expression: 

gx  =  8  /i(4>) 


where  f^)  = 


sin<j> 


(H-21) 


v/sin2^  +  /t2cos24» 


5.  Potential  Flow  at  Ellipse  Surface, 

When  considering  forced  convection  due  to  the  flow  of  vapor  over  the 
condensate  film,  the  velocity  of  the  vapor  influences  the  condensate  film  thickness 
through  the  vapor  shear  at  the  film/vapor  interface.  The  vapor  velocity  at  the  interface 
will  be  dependent  on  the  vapor  velocity  outside  the  vapor  boundary- layer.  This  velocity 
is  determined  from  potential  flow  theory.  Assuming  that  the  film  thickness,  6,  and  the 
vapor  boundary-layer  thickness,  A,  are  much  smaller  than  the  radial  distance,  r,  of  the 
elliptical  surface  (and  therefore  may  be  neglected),  the  potential  flow,  U^,  about  an 
elliptical  surface  with  major  axis  aligned  with  the  vapor  free  stream  velocity,  U„  is 
given  by 

U.  =  V-  -<1+k)  -  •  (H-22) 

\Jl  *  A^cotfy 

Details  of  the  derivation  of  this  expression  are  provided  in  Appendix  (A). 

6.  Ellipse  with  Major  Axis  Perpendicular  to  Gravity 

For  the  case  where  the  ellipse  major  axis  is  perpendicular  to  the  vertical  as 
shown  in  Figure  II-3,the  previously  determined  equation  for  streamwise  component  of 
gravity  remains  the  same.  For  0  <  k  <  1,  the  major  axis  is  aligned  with  gravity.  For  1 
<  k  <  <*>,  the  major  axis  is  perpendicular  to  gravity. 
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Figure  II - 3.  Geometry  for  Condensation  on  Horizontal  Elliptical 
Cylinder  with  Major  Axis  Perpendicular  to  Gravity. 


B.  FREE  CONVECTION  CONDENSATION  ON  AN  ELLIPTICAL  CYLINDER 
(NUSSELT  [5]  TYPE  ANALYSIS) 

Consider  a  quiescent,  saturated  vapor  at  a  temperature  Tsal,  condensing  on  a 
horizontal  elliptical  cylinder  of  eccentricity,  k,  and  semi- major  axis,  a,  aligned  with  the 
direction  of  gravity  as  in  Figure  II- 1.  The  same  simplifying  assumptions  are  made  as 
proposed  by  Nusselt  in  his  analysis  of  condensation  on  a  horizontal  circular  cylinder. 
These  assumptions  are: 

(1)  The  vapor  is  pure  and  quiescent. 

(2)  The  tube  wall  temperature  is  uniform  and  constant. 

(3)  The  thermophysical  properties  of  the  fluid  are  constant  and  evaluated  at  some 
given  reference  temperature. 

(4)  The  temperature  at  the  film/vapor  interface  is  Tsat. 

(5)  The  condensate  film  thickness  is  small  compared  to  the  radius  of  the  ellipse. 

(6)  Heat  transfer  in  the  condensate  film  is  one  dimensional  in  the  radial  direction 
providing  a  linear  temperature  distribution  across  the  film  (based  on  5«R). 

(7)  The  only  forces  acting  on  the  condensate  film  are  viscous  and  gravity  forces. 

(8)  Flow  in  the  liquid  film  is  laminar  with  no  waviness  and  no  vapor  boundary 
layer  separation. 
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A  fluid  element  ot  unit  width  is  analyzed  in  detail  and  shown  in  Figure  II-4. 


Convection. 


Conservation  of  mass  for  the  fluid  element  requires  that, 

[pi/0  u  dy  dx  -  [p,/\ dy  +  ~  jpj^u  dx  *  m  dx  -  0  (11-23) 

where  m  is  the  local  condensation  mass  flux  rate. 

Defining  the  mean  film  velocity  as 

“» *  ifou *  ’  (II-24) 

Equation  (II - 23)  simplifies  to 

«  =  (n-25) 

Conservation  of  momentum  for  the  film  element  is  a  balance  of  viscous  and  body 
forces  which  reduces  to 

1,^-f  *  9,8  /,(♦)»  0  .  (II -26) 

dy1 
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(II  -  27) 


Here,  an  assumption  is  made  that  the  fluid  is  Newtonian,  i.e., 

du 


’It 


dy 


Conservation  of  energy  is  a  balance  of  latent  heat  from  condensation  and  heat 
conduction  through  the  condensate  film  which  reduces  to 


hfjn  dx  =  dx 


m  = 


A» 


(II -28) 


V 


The  boundary  conditions  for  the  momentum  equation  are 


“,-o  =  0 


*8  = 


( du\  „ 
—  =0 

W,-* 


(11-29) 


Integrating  the  momentum  equation  subject  to  the  boundary  conditions  results  in 
a  local  film  velocity  of 

<n-30) 


and  a  mean  film  velocity  of 


(11-31) 


Combining  the  continuity  equation,  energy  equation  and  expression  for  film  mean 


velocity  (Equations  (II-25),(II-28)  and  (II- 31)),  results  in 


A,,  AT 


d_ 

P?  g 

dx 

(II -32) 


Using  the  relationship  between  dx  and  d$  from  Equations  (II-14)and  (II-15)and 
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carrying  out  the  differentiation  results  in 


AT  i  =  pf_g 

hh  5  3Tb 


fa3  X«t»)  -  3  /,«*>)  62£ 

m  i  d$ 


(II-33) 


A  dimensionless  parameter,  Zt,  is  defined  by 

‘ 


(II -34) 


where  Bj  is  given  by 


k,  ii;  A T  Dt 

P?  hf*8 


(II  -  35) 


Note  that  the  equivalent  parameter  used  by  Nusselt  [3]  multiplies  Equation  (II - 35)  by 
a  factor  of  three  and  uses  radius  in  place  of  effective  diameter.  Substituting  Equations 
(11-34)  and  (II-35)into  (11-33)  results  in  a  first  order  ordinary  differential  equation  for 

Zv 

dZ\  +  4/fo»  z  _  4  x(4>)  (11-36) 

d*  3/j(<|>)  1  DJW 


This  equation  can  be  solved  exactly  by  using  an  integrating  factor  which  would 
require  numerical  integration  or  can  be  solved  numerically  using  a  forward  stepping  or 
propagation  technique.  The  initial  condition  is  determined  based  on  the  symmetry  of 
the  problem  which  requires  that  (dZj/d^-o  =  0.  Solving  Equation  (II- 36)  at  4>  =  0 
results  in  the  initial  condition 


=  3  X(0)  =  3  a*2 


(11-37) 
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The  non-dimensional  film  thickness,  6*,  and  film  thickness,  6,  are  determined 
respectively  by 


8*(4>)  =  Zj(<t>),/4  and 

*(♦)  =  <  z,(*)  fi,  • 


(11-38) 


The  local  heat-transfer  coefficient,  a,  and  mean  heat-transfer  coefficient  are  given 
respectively  by 


a 


a(4>)  = 


m 


and 


— /“«(♦)  X(4>)  , 

7tZ)/0 


(II -39) 


and  the  local  Nusselt  number  and  mean  Nusselt  number  by 


Afo«t>)  = 


«(4>)X>« 


and 


m 


ip_L 

8(d>) 


X(4>)  d* 


(11-40) 


For  the  case  where  k  =  1  (circular  tube);  f  j(<J>)=sin<f>,  f  j  '(4*)=cos<f>,  %(<f>)=a  and  De  = 

2a  and  Equation  (II -36)  reduces  to 

<Ei  +  4cos4»  z  =  JL  f  (11-41) 

<?<(>  3  sin4>  1  sin<j> 

which  is  the  same  as  that  found  by  Nusselt  [5]  except  for  a  difference  in  the  definition 
of  Br 

The  above  development  can  also  be  applied  to  a  horizontal  elliptical  tube  whose 
major  axis  is  perpendicular  to  the  direction  of  gravity  as  shown  in  Figure  II -3.  In  this 
case,  b  is  greater  than  a,  i.e.  k  >  1. 
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C.  FORCED  CONVECTION  CONDENSATION  ON  AN  ELLIPTICAL  CYLINDER 
(SHEKRILADZE-GOMELAURI  [21]  TYPE  ANALYSIS) 

Consider  a  saturated  vapor  at  temperature  Tsat,  flowing  downward  over  a 
hori7ontal  elliptical  cylinder  with  free  stream  vapor  velocity,  U„.  The  elliptical 
cylinder  has  eccentricity,  k,  and  semi-major  axis,  a,  parallel  with  the  direction  of 
gravity  and  vapor  flow.  The  same  assumptions  as  for  the  case  of  free  convection  are 
used  here  with  the  exceptions: 

(1)  The  force  of  gravity  is  neglected. 

(3)  U,»uy=5. 

(4)  The  interfacial  shear  stress  is  approximated  by  an  asymptotic  expression 
assuming  an  infinite  condensation  rate, 

=  mU *  (II -42) 


Conservation  of  mass  and  energy  are  the  same  as  previously  derived  (Equations 
(11-25)  and  (11-28)).  The  only  forces  acting  on  the  fluid  element  are  viscous  forces. 
Thus,  conservation  of  momentum  reduces  to 

d2u 


with  the  boundary  conditions, 


’ll  —2  =  0 
<ty 


Uy. 0  =  0 


l  du\  ,, 

1i  j  - 

\  dy  )y.t 


(11-43) 


(II -44) 


Integrating  the  momentum  equation  subject  to  the  boundary  conditions  results  in  a 
local  velocity  of 


«(y)  = 


m  U. 


’ll 


(II-  ') 
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and  a  mean  film  velocity 


u_  = 


Hi  2 


(II -46) 


Using  the  expression  for  potential  flow  over  an  elliptical  cylinder  (Equation  (II-22))and 
defining  a  velocity  function,  f2(4>),  such  that 


m)  = 


1+* 


/l  +  ifc2cot2<j) 


(11-47) 


the  mean  film  velocity  can  be  expressed  as 

m  t/_ 


u_  = 


*1/  2 


(11-48) 


Combining  the  continuity  equation,  energy  equation  and  expression  for  mean  film 

velocity  (Equations  (II-25),(II-28)  and  (II- 48))  results  in 

kt  A7~  j.  _  d  \ Pi  2,  A7-  Um  (II- 49) 

\  8  *1  Hi**  2  J  • 

Using  the  relationship  between  dx  and  df  from  Equations  (II -14)  and  (11-15),  and 
carrying  out  the  differentiation  results  in 


k,  A  T  i  p,  X,  A  T  U.  i 


* ft  8 


2t1  iK 


Defining  a  dimensionless  parameter,  Z2,  as 

z*  - 


6^ 

*2 


(H-51) 


where 


*2  = 


n  iD. 

p,  v. 


ai-52) 
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Equation  (11-50)  becomes 


Z,  -  4  (II -53) 

d<t>  W)  ^  D./20 t>) 


This  problem  is  a  first  order  ordinary  differential  equation  which  can  be  solved  using 
a  forward  stepping  numerical  integration  technique.  The  initial  condition  for  the 
problem  is  based  on  symmetry  which  requires  that  (dZ2/d<|>)^_0  =  0.  Solving  Equation 
(II -53)  at  <j>  =  0  results  in  the  initial  condition 


=  2  X(0)  .  2a k2 
Djm  °e(1+*> 


(11-54) 


The  non-dimensional  film  thickness,  6  ,  and  film  thickness,  6,  are  determined 
respectively  by 


**  = 

8  =  • 


(11-55) 


The  heat- transfer  coefficient  for  forced  convection  condensation  can  be  considered  in 
a  different  form.  Combining  Equations  (II-40),(II-51),  (11-52)  and  (II -55)  results  in 


Nu  Re  1/2  =  — - —  ,  and 
Nu  Re  1/2  =  —/* — ~ —  X(4>) 


(II -56) 


For  the  case  where  k  =  1  (circular  tube),  /2(<f>)  =  2sin 4>,  dx  =  nd<J>  and  Equation  (II- 
49)  reduces  to 


6 


1  d  P<  U. 
ad<J>[  q, 


sin«j>  6 


(II -57) 


which  is  the  same  as  the  Shekriladze-Gomelauri  [21 J  governing  equation  for  a  circular 
cylinder  with  no  body  forces. 
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D.  MIXED  CONVECTION  CONDENSATION  ON  AN  ELLIPTICAL  CYLINDER 

(SHEKRILADZE-GOMELAURI  [21]  TYPE  ANALYSIS) 

Consider  a  saturated  vapor  at  temperature  Tsat,  flowing  downward  over  a 
horizontal  elliptical  cylinder  with  free  stream  vapor  velocity,  U„.  The  elliptical 
cylinder  has  eccentricity,  k.  and  semi-major  axis,  a,  which  is  parallel  with  the  direction 
of  gravity  and  vapor  flow.  The  forces  acting  on  a  film  eltment  are  viscous  and  body 
forces.  The  interfacial  shear  is  given  by  the  Shekriladze-Gomelauri  [21]  model  for 
infinite  condensation  rate  as  in  the  case  of  the  forced  convection  model  (Equation  (II- 
42)). 


Conservation  of  mass  and  energy  are  as  previously  derived.  Conservation  of 
momentum  for  the  condensate  film  reduces  to 

(II -58) 


’IiT 1  +  P  !*/.<♦>  =° 
<fy 


with  boundary  conditions, 


“,.o  -  0 

Mu)  =  m  U* 

V  dy  Jj„4  q( 


(11-59) 


Integrating  the  momentum  equation,  subject  to  the  boundary  conditions,  results 
in  the  local  film  velocity 


Pi  8  (  v2^  m 

«(y)  = - Uy-^-  + - /2(<t>)  y 

*1/1  2  /  q, 


(II -60) 


and  mean  film  velocity 


Pi  8  ft2  m  U_  a 

=  ~  -  — -m)  f 

3  q}  2 


(II -61) 


where  the  expression  for  potential  flow  over  an  ellipse  and  its  associated  velocity 
function  (Equations  (II -22)  and  (11-47))  have  been  used. 
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Combining  the  continuity  equation,  energy  equation  and  expression  for  mean  film 


velocity  (Equations  (11-25),  (11-28)  and  (II-61))  results  in 

*1  Ar  ±  =  A pfg  8*  +  P,  X,  AT  U, 
/tA  6  dxjq,  1  3  rt(AA 


(11-62) 


Using  the  relationship  between  dx  and  dip,  and  algebraically  manipulating  Equation  (II- 
62)  results  in 


_1_ 

6* 


2  xft) 


2 *7,(0)  -y 


+  m)  y 


(11-63) 


where  F  is  a  dimensionless  parameter  relating  free  and  forced  convection  and  6*  is  a 
dimensionless  film  thickness  given  respectively  by 


F  = 


-  ^  D«  hh  8 


U'\,  A  T 


y  = a 


n  r 


n, 


(II -64) 


Differentiating  Equation  (11-63)  results  in 


_1_ 

5* 


D. 


2  x(4>) 


2F /,(<!>)  6*1  ♦  2F/fa>) 


(11-65) 


Regrouping  the  terms  shows  that  the  equation  is  a  first  order  ordinary  differential 
equation: 


{2F/,(4>)  A  *  m  5*}^ 


-  2  X(4» 

D, 


(II -66) 


The  initial  condition  for  this  problem  is  again  based  on  symmetry,  which  requires  that 
d6*/d<|>=0  at  <J>=0.  Inserting  this  condition  into  Equation  (II -66)  results  in  a  fourth  oru  r 
expression  for  5* 

I  F/fc 0)  y*  *  f2(0)  yl  -  2  =  0  ,  (11-67, 


38 


whose  root  of  interest  is  given  by 


«'( 0)  = 


+  4  <2/3  F/f(0)}  {2 


2  (2/3  Ff[(0)) 


3  k 

(Uk\\ 

16  F  a 

1+* 

4  F 

\|l  *  i 

3D. 

k 

/J 

(11-68) 


equation  (11-66)  can  be  solved  numerically  to  obtain  the  film  thickness.  The  other  heat 
transfer  parameters  are  determined  as  before  (Equations  (II- 39)  and  (II- 56)). 

For  the  case  where  F=0  and  k= 1  (forced  convection  over  a  circular  tube),  Equation 
(11-63)  reduces  to  Equation  (II-57)  as  expected.  For  the  case  where  F  -  “and  k  =1  (free 
convection  over  a  circular  tube),  Equation  (11-63)  becomes 


5* 


B  d_ 

2x(4>) 


j!  Fft( 40  a 


(11-69) 


Upon  substitution  for  the  parameters  F  and  3*  (Equation  (11-64)),  Equation  (11-69) 
becomes 


l.iAiitl  id***! 

6  a  X,  AT  r\,  d${  3 


(II -70) 


which  is  the  same  as  Nusselt’s  governing  equation  [5]  for  laminar  film  condensation  on 
a  horizontal  circular  cylinder. 


E.  MIXED  CONVECTION  CONDENSATION  ON  AN  ELLIPTICAL  CYLINDER 
WITH  SURFACE  TENSION  AND  PRESSURE  GRADIENT  EFFECTS 
Consider  a  pure  saturated  vapor  at  temperature  T^,,  flowing  downward  over  a 
horizontal  elliptical  cylinder  with  free  stream  vapor  velocity,  U„.  The  elliptical 
cylinder  has  eccentricity,  k,  and  semi -major  axis,  a,  parallel  with  the  direction  of 
gravity  and  vapor  flow.  In  addition  to  the  assumptions  used  for  mixed  convection  in 
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the  last  section,  pressure  effects  in  the  condensate  film  are  considered  in  the 
development  of  the  governing  equations.  The  pressure  on  a  condensate  film  element 
will  be  a  combination  of  that  impressed  upon  the  element  by  the  potential  flow  of  the 
vapor  as  well  as  that  due  to  surface  tension.  A  fluid  element  is  shown  in  detail  in 
Figure  II - 5. 


Conservation  of  mass  and  energy  are  as  previously  described.  Conservation  of 
momentum  for  the  film  element  is  a  balance  of  viscous,  pressure  and  body  forces  which 
reduces  to 

*,—  +  (> ,*/,(♦> (II-7D 
dyi  dx 


with  boundary  conditions 


“>-o 


=  0 


fdu)  =  "  V* 


(11-72) 
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As  previously  stated,  the  film  pressure  is  a  combination  of  the  pressure  due  to 
potential  flow,  p^,  and  the  pressure  due  to  surface  tension,  p5.  Thus  the  total  pressure 


gradient  is  given  by 

dp  _  dp  o 

dx  dx  dx 

The  pressure  gradient  due  to  potential  flow  is  given  by 

dp±  v  dU*  -pvt/«,  dU t 
dx  '  "Pv  ♦  dx  X(<t») 


(11-73) 


(11-74) 


Substituting  the  expressions  for  potential  flow  and  x(4>)  (Equation  (11-22)  and  (11-15)), 


into  Equation  (11-74)  results  in 


dx 


Pv  ^ 
la 


/,(♦) 


(11-75) 


where 


cos2^  +  fe2sm2<t> 
it2  +  ^(Jt2-1)2  sin^fj* 

The  pressure  due  to  surface  tension  is  given  by 

P  =  —  . 

’  *♦) 

It  is  assumed  that  8«R.  Thus,  the  pressure  gradient  is  given  by 

dP,  _  __o_  dR  q  dR 

dx  r}  dx  s*  x(4>)  d$ 


m  - 


(1+fc)2  *2  sin24> 
sin24>  +  fc2cos2<|»  ] ' 


(II -76) 


(II -77) 


(II -78) 


Substituting  Equations  (II -9)  and  (11-15)  into  Equation  (II- 78)  results  in 


dP, 

dx 


(II -79) 
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where 
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dimensionless  parameters  defined  by: 


p  _  py  hft  ** i 

P, 

P'S* 


(11-85) 


P  is  related  to  the  pressure  gradient  and  Bo,  the  Bond  number,  relates  the  inertia  effects 
to  surface  tension  effects.  Differentiating  Equation  (11-84)  results  in  the  following  first 
order  ordinary  differential  equation: 

=  2  F/,( 40  +  ^  F/3(4>)  -  3  (-*)  /4(<t>)  6*J 

De  ( [  a  a  )  Bo 

j»*  f  .  D  ,  (11-86) 

+  /2(<t>)  »•}-£-  +  2  F/fa)  +  —  Pfo) 

09  (  <2 

In  the  previous  development,  film  separation  and  instability  were  not  applicable 
since  no  effects  were  considered  which  could  retard  the  film  flow  and  (dfi*/d<J>)-<»only 
at  However,  in  the  present  analysis,  the  decreasing  potential  flow  and  increasing 

radius  of  curvature  of  the  elliptical  cylinder  causes  the  pressure  gradient  to  have  a 
retarding  effect  on  the  condensate  film  over  the  back  half  of  the  tube.  At  the  point  of 
condensate  film  separation  from  the  tube  wall,  the  film  shear  stress  at  the  wall  is 
negative  which  implies  flow  reversal  (du/dy  s0aty  =  0).  Differentiating  Equation  (II- 
81)  and  solving  for  the  film  separation  criteria  (du/dy  =  0  at  y  =  0)  results  in  the 


following  condition  for  separation  to  occur: 

2F /j(4>)  *  -!  P/3(<t>)  +  3  (^)  -£-  /<($)  6*’  +  2/2(<|>)  6*  s  0 
a  \  a  Bo 


(II -87) 
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Condensate  film  instability  occurs  when  the  rate  of  change  of  film  thickness  is  infinite. 
From  Equation  (II -86),  this  condition  occurs  when 

D  (d\2  f 

2F /j(4>)  -  —  Pm)  +  3  -J  -jL  /4(<t») 
a  \  a  )  Bo 

From  a  comparison  of  Equations  (II -87)  and  (11-88),  it  is  evident  that  the  condition  of 
film  instability  will  occur  first  (has  a  smaller  magnitude  than  the  film  separation 
criteria  for  all  conditions)  and  is  therefore  the  limiting  condition.  Rose  [28]  suggests 
that  this  instability  may  manifest  itself  in  wavy  or  turbulent  flow,  or  may  signify  the 
detachment  of  the  liquid  film  from  the  tube  surface.  For  the  case  where  1/Bo  =  0 
(surface  tension  neglected),  and  using  Equations  (11-21), (II-47)and  (11-76)  for /j(<J>), /,(<(>) 
and  /3( <J>)  in  Equation  (11-88),  it  can  be  seen  that  the  condition  for  which  solutions  can 
be  determined  over  the  entire  tube  surface  occurs  when 

P  <  — -)2  F  (II -89) 

D\Uk) 

For  a  circular  tube  (fc= 1),  Equation  (11-89)  reduces  to  P<F/ 8,  which  was  that  found  by 
Rose  [28].  If  the  opposite  of  Equation  (11-89)  is  true,  then  the  pressure  gradient  effect 
is  dominant  enough  to  significantly  retard  the  flow  over  the  back  half  of  the  tube  and 
result  in  a  rapid  thickening  of  the  condensate  film.  Note  that  the  condition  for 
condensate  film  instability  is  a  function  of  vapor  velocity  ( F )  since  the  pressure 
gradient  due  to  potential  flow  is  velocity  dependent. 

For  the  case  where  P  =  0,  and  using  Equation  (II -  88),  the  condition  which  allows 
solution  of  the  momentum  equation  over  the  entire  tube  is  now  given  by 

JL  <  i  f-g-f  k*  (II  -9i  ' 

Bo  3  (dJ  l -*2 


y  +/2(4>)  y  =o 


(II -88) 
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Note  here  that  the  condition  for  film  instability  is  only  a  function  of  geometry. 
Therefore,  for  small  values  of  eccentricity  (approaching  a  vertical  flat  plate),  the 
retarding  effect  of  surface  tension  on  the  condensate  film  over  the  back  half  of  the 
tube  is  significant  and  can  result  in  a  rapid  thickening  of  the  condensate  film  at  <J>C. 

The  initial  condition  for  this  problem  is  again  based  on  symmetry,  which  requires 
that  (dfi’/d<t»)^_0  =  0.  Combining  this  condition  with  Equation  (11-86)  results  in  a  fourth 
order  expression  for  6  : 


2  F/fc 0)  +  ^  P/3(0)  *  3  (^)  /4(0) 

2a  a  )  Bo 


~  +  f2(0)  6* 


2  *(°)  =  o  (H-91) 


whose  root  of  interest  is  determined  by 

2F/l(0)  +  ^  P/3(0)  *  3  f— |  ~  fA( 0) 
[  a  \  a  Bo 

B  =/2(0) 


A-h 

3 


C  = 


_2X(0) 


(11-92) 


«*(0)  = 


y/B2-4  AC  -  B 
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Equation  (II -86)  can  be  solved  numerically  to  obtain  the  condensate  film  thickness.  The 
other  heat  transfer  parameters  are  determined  as  before  (Equations  (11-39)  and  (II - 56)). 

It  can  be  shown  that  this  formulation  reduces  to  the  models  developed  for  circular 
tubes.  Consider  the  case  where  k= 1  and  l/Bo=0  (circular  tube  with  no  surface  tension); 
Equation  (11-84)  reduces  to: 


_1_ 

6* 


ft  * 

2  F  sin<t>  +  8  P  sin24>  ]  —  +  2  5'  sin4> 


(II  -  93) 
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which  is  identical  to  that  developed  by  Rose  [28]  where 


6’  =  6 


( «-  0.) 


(2d  Tl, 


v/2 


y 


feat 


(II -94) 


For  the  case  where  Jfc=l,  P=0  and  l/f?o=0  (circular  tube,  no  surface  tension  and  no 
pressure  gradient),  Equation  (11-84)  reduces  to  the  Shekriladze-Gomelauri  [21]  model. 


F.  MIXED  CONVECTION  CONDENSATION  ON  AN  ELLIPTICAL  CYLINDER 
WITH  VAPOR  BOUNDARY  LAYER  SEPARATION  (FUJII  ET  AL.  [22J  TYPE 
ANALYSIS) 

The  mixed  convection  model  developed  by  Shekriladze  and  Gomelauri  [21]  has 
been  used  because  of  its  simplicity  anH  ease  of  solution.  However,  since  the  interfacial 
shear  is  approximated  by  an  asymptotic  expression  based  on  an  infinite  condensation 
rate  and  uses  potential  flow  outside  the  vapor  boundary -layer  (which  is  always 
positive),  vapor  boundary -layer  separation  does  not  occur.  As  described  in  Chapter  I, 
Fujii  et  al.  [22]  modified  the  interfacial  shear  stress  expression  by  simultaneously 
solving  the  boundary -layer  equations  for  the  condensate  and  vapor  ensuring  compatibil¬ 
ity  at  the  condensate/ vapor  interface.  This  technique  allows  for  a  more  precise 
description  of  the  condensation  problem  under  forced  convection  conditions  and 
includes  vapor  boundary -layer  separation. 

Consider  a  pure  saturated  vapor  at  temperature  Tsat,  flowing  downward  over  a 
horizontal  elliptical  cylinder  with  free  stream  velocity,  U..  The  elliptical  cylinder  is 
oriented  such  that  its  major  axis  is  aligned  with  the  direction  of  gravity  and  vap  flow. 
Due  to  the  complexities  of  solving  the  two-phase  boundary- layer  equation,  it  is  as  led 
that  the  only  forces  acting  on  the  condensate  film  element  are  gravity  and  v.  is 
forces.  Surface  tension  and  pressure  gradient  are  neglected  which  is  reasonable  s  e 
they  were  found  to  have  negligible  influence  under  most  conditions  (discusset  n 
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Chapter  IV).  Other  assumptions  are  the  same  as  in  Section  II. E.  The  vapor  is  assumed 
to  be  in  thermal  equilibrium  such  that  the  energy  equation  for  the  vapor  need  not  be 
considered  in  the  analysis.  The  only  forces  acting  on  an  element  in  the  vapor  boundary- 
layer  (see  Figure  II - 6)  are  inertia,  viscosity  and  pressure  gradient,  which  is  impressed 
upon  the  element  by  potential  flow  outside  the  vapor  boundary- layer.  The  tangential 


Figure  II-6.  Condensate  Film  and  Vapor  Elements  for  Mixed  Convection  using 
Fujii  et  al.  [22]  Type  Model. 


velocity  at  the  vapor/condensate  interface  is  much  smaller  than  the  potential  velocity 
outside  the  vapor  boundary- layer  and  can  therefore  be  assumed  to  be  negligible.  Based 
on  the  above  assumptions,  conservation  of  mass  and  momentum  for  the  vapor  boundary- 
layer  are  given  by 

a//  ai/ 

( continuity )  (a) 


dx 


aUL.vSL.n 


dU  dV  - 

—  +  —  =o 


dy 

dU, 


»  ,  n„  &U 


(11-95) 
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dy 


dx 


P  y  dy 


( momentum )  (b) 
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For  the  condensate  film,  conservation  of  mass,  momentum  and  energy  are  given  by 


—  +  g  /j(4>)  =  0  ( momentum )  (a) 

P,  dy2 

M  =  “*  ^  =  - 1 — * -  ( ?”*rgy  and  mass )  (b) 


(II -96) 


h A* 


where  conservation  of  mass  and  energy  have  been  combined  as  in  previous  sections. 
Boundary  conditions  for  this  system  of  equations  are  given  by 


«y-o  =  Vo  =  0 


U  =  u  (—) 
♦  ’  (  dy  )y 


(II -97) 


=  0  . 


6*6. 


where  A  is  the  vapor  boundary -layer  thickness.  The  compatibility  conditions  at  the 
interface  require  that 

*1  (Tam~Twati} 


(11-98) 


"PvK,-a  = 


hft* 


where  the  last  compatibility  relation  is  the  condensation  mass  flux  rate  (i.e.  for  a  given 
surface  area,  the  amount  of  vapor  condensed  must  equal  the  increase  in  the  mass  of  the 
film). 

The  technique  for  solving  the  vapor  boundary-layer  momentum  equation  involves 
an  approximate  integral  technique  developed  by  Truckenbrodt  [23J.  The  initial  problem 
development  follows  the  method  of  Pohlhausen  [36],  The  continuity  and  momentum 
equations  (Equations  (II-95a)  and  (II- 95b)),  are  integrated  over  the  thickness  of  the 
boundary  layer  from  y  =  6  to  y  =  d+A  and  then  combined  to  eliminate  the  normal 
component  of  vapor  velocity,  V,  resulting  in  a  momentum  integral  equation 

(11-99) 


±  A,}  *  ^  „  _lL 

ul  *'  *  "  V,  dx  vt  p>  V{ 
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where  Ajand  A,  are  the  vapor  boundary  layer  displacement  and  momentum  thicknesses 
given  respectively  by 


A, 


and  A2 


U_ 
U A 


(II -100) 


V5  is  the  normal  component  of  vapor  velocity  at  the  vapor/condensate  interface  due  to 
the  condensation  process  itself  and  is  sometimes  called  the  vapor  boundary -layer  suction 
velocity.  Differentiating  the  A2  term  in  Equation  (11-99)  and  rearranging  the  equation 
results  in: 


<*2  ,  dU*  2  .  A1  )  _  V>  .  ^ 
djc  dx  A,  J  Py 

With  some  further  algebraic  manipulation  this  yields 


d  < 

K  pvj 

2 

A2 

a2  Pv  du ;  Kj  a2  pv 

dx] 

Hv  j 

tlv  ^ 

r 

*iv  *  nv 

(II- 101) 


(II- 102) 


To  solve  this  expression,  the  vapor  velocity  distribution  across  the  boundary- layer 
thickness  must  be  determined.  This  distribution  is  a  function  of  a  pressure  gradient 
parameter,  k,  and  a  suction  parameter,  Kj,  given  by  Tuckenbrodt  [23]  as 

V%  *2  Pv 


*1  Pv  , 

x  = - -  and 

1\y  dX 


Kj  =  - 


(II  103) 


Truckenbrodt  [23]  defined  a  shear  function  f(K,Kj) 


uuu 


/fo^)  = 


T»  A2 
*1,  U+ 


and  H(k,  k,) 


(11-104) 


and  a  dimensionless  function,  Z,  given  by 


Z  = 


(II -105) 
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Equation  (II -102)  becomes 
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The  initial  condition  for  Equation  (11-112)  conies  from  symmetry  and  is  given  by 

(—)  =0  .  (11-113) 

Combining  Equations  (II-112)and  (11-113)  results  in 

6  [ ^*\  L/ZfjO)  f  *  0.882  G  —(0)  s[ZQ)  ~  0  441  =  0  . 

X(0)  l  d< t>  y 


(11-114) 


This  is  a  quadratic  equation  from  which  Z( 0)  is  determined  to  be, 
Z( 0)  = 


-0.0735  G  x(0) 

i 

l<*4>  V> 

6*(0) 

0.0735  G 


X(0) 


T  .  0.0735  m 


(II- 115) 


Vo  6  *(0)  J  De  \  d*  Jtm0 

The  interfacial  shear  is  determined  from  Equation  (II- 104).  Truckenbrodt  [23] 
provides  a  simple  approximation  for  the  shear  function,  /(k.Kj),  given  by 

/(K,K,)  =  3.22  fa  (Ka+K)  (H-  H6) 

where 

k  =  0.0682  +  0.174  k, 


k  and  Kj  are  modified  by  combining  Equations  (II-103)and  (II- 105): 

Dt  dU*  ,  G  r- 

K  =  _i - ♦  Z  and  k.  =  —  yfz  . 

X(4>)  d< J)  1  5* 

Dimensionless  interfacial  shear  is  defined  by: 

_  2  r4 
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(11-117) 


(11-118) 


(II- 119) 


An  expression  for  dimensionless  shear  can  then  be  obtained  by  combining  Equations  (II- 
104),  (II -105)  and  (II -116): 


=  6.44-*  . 

y/Z 


(11-120) 
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The  governing  equations  for  the  condensate  film  are  solved  in  the  usual  manner. 
The  momentum  equation  (Equation  (II -96a)),  is  integrated  to  solve  for  condensate  film 


velocity, 


u 


The  mean  velocity  is  determined  from  Equation  (11-24): 


u 


M 


Pi  8 


(11-121) 


(11-122) 


Substituting  Equation  (11-122)  into  the  energy/continuity  equation  (Equation  (II-96b)) 
and  multiplying  both  sides  of  the  equation  by 

n'  D'  hf* _  (II  - 1 23) 

Pi  Ul  X, 


results  in 
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x(4»)  d* K 


+  I±7 
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(II- 124) 


Carrying  out  the  differentiation  of  Equation  (11-124)  and  rearranging  results  in  the 


ordinary  differential  equation, 


dr 

d«(» 


X(4>)  F  dfx  y*  .  1  1  dxh  8.3 

Dt  3  d<t>  4  G  d<p 


(II- 125) 


which,  due  to  the  symmetry  of  the  problem,  has  the  initial  condition, 


V  d*  l4+ 


(II -126) 
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Combining  Equations  (11-125)  and  (II- 126)  results  in 


X(0)  £ 
Dt  3 


d*‘(0)  -  -  - 
„  4  G 


6  *3(0)  =  0 
>=o 


(II- 127) 


from  which  6  (0)  can  be  determined  by  numerical  methods. 

For  the  case  of  k  =  1  (circular  tube),  Equation  (II- 125)  reduces  to  the  analysis  of 
Fujii  et  al.  [22]  as  corrected  by  Lee  and  Rose  [25].  The  differential  equation  for  the 
vapor  boundary- layer  (Equation  (11-112))  and  for  the  condensate  film  (Equation  (II- 
125))  are  solved  simultaneously  using  the  compatibility  relations  of  Equations  (II- 117), 
(II-118)and  (11-120).  At  the  point  of  vapor  boundary-layer  separation,  the  interfacial 
shear  stress  becomes  negative.  Downstream  of  this  point,  it  is  assumed  that  the  shear 
stress  at  the  vapor/condensate  interface  is  negligible  and  a  simple  Nusselt  type  analysis 
is  used  for  the  remainder  of  the  elliptical  surface.  The  heat  transfer  parameters  are 
determined  as  before  (Equations  (II -39)  and  (11-56)). 
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III.  NUMERICAL  METHOD 

In  the  theoretical  development,  the  governing  equations  for  each  model  were 
either  reduced  to  a  single  first  order  ordinary  differential  equation  (ODE)  or  to  a 
simultaneous  solution  of  a  system  of  first  order  ODEs  of  the  form 

&  =  M*J)  ■  (IM-l) 

dx 

Initial  conditions  were  determined  from  the  symmetry  of  the  problem  which  required 
that  both  the  slopes  of  the  condensate  film  and  vapor  boundary- layer  wen.  zero  at  <j>  = 
0.  These  systems  of  equations  are  solved  using  one  of  the  forward  stepping  numerical 
methods. 

A.  SOLUTION  TECHNIQUE  FOR  THE  NUSSELT  [S]  AND  SHEKRILADZE- 

GOMELAURI  [21]  ANALYSIS  METHODS 

Each  analysis  method  resulted  in  a  single  first  order  ODE  of  the  form  given  by 
Equation  (III- 1).  The  solution  of  this  equation  provides  the  condensate  film  thickness 
along  the  elliptical  tube  surface.  Many  numerical  methods  are  available  to  solve  this 
type  of  problem,  each  with  its  own  advantages  and  disadvantages.  For  this  study,  a 
fifth  order  Adams  method  predictor-corrector  algorithm  was  developed  from  Crandall 
[39].  The  Adams  methods  are  multi-step  techniques  which  use  finite  difference  type 
operators  incorporating  previously  determined  points  (hence  the  term  multi-step)  to 
determine  the  value  corresponding  to  the  next  step.  These  methods  provide  the  same 
accuracy  and  are  more  efficient  than  the  Runge-Kutta  (single  step)  methods  but  re 
restricted  to  a  fixed  step  size.  The  Adams  method  employed  in  this  study  uses  an 
Adams-Bashforth  explicit  method  [39]  to  predict  the  value  of  the  dependent  variable  at 
the  next  step  and  then  an  Adams- Moulton  implicit  method  [39]  to  correct  this  predicted 
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value.  This  method  is  not  self  starting  and  requires  determination  of  the  first  four 
points.  To  start-up  the  algorithm,  lower  order  Adams  methods  were  used  to  obtain  these 
four  points.  The  recurrence  formulae  for  this  algorithm  are  given  by 

y,  =  Vi  +  h  (in- 2) 

k 

where  x  is  the  independent  variable,  y  is  the  dependent  variable,  h  is  the  step  size,  s  is 
the  step  number  for  which  y  is  to  be  determined,  k  are  the  points  determined  in 
previous  steps  and  Pk  are  the  coefficients  provided  in  Table  III  - 1  and  III-2. 


TABLE  III- 1.  RECURRENCE  COEFFICIENTS  FOR  PREDICTOR  FORMULA 


FORMULA 

Pk 

s-5  s-4  s-3  s-2  s-1 

TRUNCA¬ 

TION 

ERROR 

(1) 

1 

0(h2) 

(2) 

-1/2  3/2 

0(h3) 

(3) 

5/12  -16/12  23/12 

o(h4) 

(4) 

-9/24  37/24  -59/24  55/24 

o(h5) 

TABLE  III-2. 

RECURRENCE  COEFFICIENTS 

FOR  CORRECTOR 

FORMULA 

FORMULA 

Pk 

s-4  s-3  s-2 

s- 1 

s 

TRUNCA¬ 

TION 

ERROR 

(1) 

1/2 

1/2 

0(h3) 

(2) 

-1/12 

8/12 

5/12 

0(h4) 

(3) 

1/24  -5/24 

19/24 

9/24 

o(h5) 

(4) 

-19/720  106/720  -264/720 

646/720 

251/720 

o(h6) 

The  algorithm  used  in  the  solution  for  mixed  convection  condensation  with  surface 
tension  and  pressure  gradient  (Section  II.E)  is  provided  in  Appendix  (B). 
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For  many  applications,  this  type  of  start-up  routine  would  normally  not  provide 
sufficient  accuracy.  However,  no  difficulties  were  encountered  using  this  numerical 
procedure  for  solving  the  asymptotic  interfacial  shear  stress  models.  Solution 
convergence  was  obtained  for  aa  angle  step  size  of  0.1°.  The  difference  between  an 
angle  step  size  of  1.0°  and  0.1°  was  less  than  1%. 

B.  SOLUTION  TECHNIQUE  FOR  THE  FUJII  ET  AL.  [22]  ANALYSIS  METHOD 

Reduction  of  the  governing  equations  resulted  in  a  system  of  two  simultaneous 
ODEs,  one  for  the  momentum  thickness  of  the  vapor  boundary- layer  and  one  for  the 
condensate  film  thickness.  The  Ad_ms  method  previously  described  did  not  provide 
sufficient  accuracy  to  solve  this  two- phase  boundary  layer  problem.  In  particular,  the 
system  of  equations  exhibited  a  stiffness  problem,  i.e.,  the  numerical  method  could  not 
accurately  solve  the  vapor  boundary-layer  ODE  over  the  first  several  steps  even  with 
an  angle  step  size  of  0.001°.  To  obtain  sufficient  accuracy,  a  much  smaller  step  size 
would  be  required  which  would  significantly  reduce  the  efficiency  of  the  computer 
algorithm.  Lee  [40J  indicated  that  similar  difficulties  had  been  encountered  when  Lee 
and  Rose  [25]  solved  the  same  equations  for  a  horizontal  circular  cylinder.  Lee  and 
Rose  [25]  used  a  Runge-Kutta  method  which  permitted  variable  stcn  size.  Thus,  the  step 
size  could  be  reduced  sufficiently  to  obtain  the  required  accuracy  for  the  vapor 
boundary-layer  ODE  over  the  front  of  the  tube  and  then  increased  over  the  remainder 
of  the  tube  to  improve  the  efficiency  of  the  algorithm.  Utilizing  this  numerical, 
procedure,  they  were  able  to  solve  the  system  of  equations  over  most  of  the  range  of 
dimensionless  parameters  ( F  and  O').  However,  in  some  cases  (specifically  large  values 
of  F  and  G),  the  system  of  equations  were  still  too  stiff  for  solution. 

To  overcome  the  problem  of  stiffness,  a  numerical  problem  solver  (IVPAG)  from 
the  International  Mathematical  and  Statistical  Library  (IMSL)  [41]  was  used.  This 
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algorithm  uses  Gear’s  stiff  method  [42]  to  efficiently  solve  systems  in  which  stiffness 
is  a  problem.  However,  as  in  the  study  of  Lee  and  Rose  [25],  this  method  could  still  not 
solve  the  two-phase  boundary -layer  equations  over  the  entire  range  of  dimensionless 
parameters  considered,  but  did  cover  the  practical  ranges  of  the  parameters. 

It  was  originally  conceived  that  the  computer  algorithm  would  progress  around 
the  tube  surface  until  the  criteria  for  vapor  boundary-layer  separation  was  reached 
(dT$/d4>  <  0).  At  this  point,  the  model  would  shift  from  solution  of  the  two-phase 
boundary- layer  equations  to  the  simple  Nusselt  model.  However,  the  numerical  solver 
encountered  difficulties  in  solving  the  system  of  equations  as  it  approached  the  point 
of  vapor  boundary -layer  separation.  Therefore,  the  solution  method  was  changed  to  a 
two  step  procedure.  In  the  first  step,  the  system  of  equations  were  solved  until  the 
computer  algorithm  indicated  a  problem  had  been  reached  (usually  the  algorithm 
indicated  a  stiffness  problem  had  been  encountered).  The  interfacial  shear  stress  data 
was  analyzed  to  verify  that  this  problem  was  the  result  of  rapidly  changing  conditions 
associated  with  vapor  boundary -layer  separation.  Therefore,  to  facilitate  the  solution, 
vapor  boundary -layer  separation  was  assumed  to  have  occurred  if  the  interfacial  shear 
stress  at  the  last  point  obtainable  was  less  than  ten  percent  of  the  maximum  interfacia) 
shear  stress  observed  over  the  surface.  At  this  point,  the  problem  solver  was  within  one 
to  two  step  sizes  of  the  actual  separation  point  (for  a  step  size  of  0.1°,  the  computed 
separation  point  was  therefore  within  0.2°  of  the  actual  separation  point).  The  second 
step  involved  re-running  the  analysis  using  the  computed  separation  point  calculated 
in  the  first  step  as  the  point  at  which  the  two- phase  model  was  switched  to  the  simple 
Nusselt  model.  A  converged  solution  was  obtained  for  an  angle  step  size  of  0.1°.  The 
algorithm  for  solving  the  two- phase  model  is  provided  in  Appendix  (C). 
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IV.  RESULTS  AND  DISCUSSION 

When  comparing  solutions  for  an  elliptical  tube  with  those  of  a  circular  tube, 
equivalent  surface  areas  have  been  used.  The  effective  diameter  of  the  equivalent 
circular  cylinder  is  given  by  Equation  (II - 16).  The  dimensionless  streamwise  length  jc’, 
defined  as  the  ratio  of  the  streamwise  length,  x,  to  the  half  perimeter  length  and  given 
by 

x*  =  -ii-  ,  (IV-1) 

7t  De 

is  used  when  comparing  local  values  of  film  thickness  and  heat-transfer  coefficients. 
Its  use  enables  direct  comparison  between  elliptical  and  circular  tubes  since  it  represents 
an  equivalent  fraction  of  the  total  perimeter  from  the  top  of  the  tube  to  the 
dimensionless  point. 

A  practical  range  of  eccentricities  for  an  elliptical  tube  with  major  axis  aligned 
vertically  is  0.3  <  k  <  0.6  .  The  lower  limit  of  0.3  is  based  on  discussions  with  a  tube 
manufacturer  (Reference  [43])  in  which  the  elliptical  tube  is  formed  by  pressing  a 
circular  tube  with  a  roller  assembly  to  obtain  the  proper  major  and  minor  axis 
dimensions.  Roller  contact  points  would  be  evenly  spaced  to  maintain  the  elliptical 
curvature  and  to  prevent  an  overly  flat  surface.  The  upper  limit  is  chosen  as  the 
eccentricity  in  which  noticeable  effects  were  observed  in  the  heat  transfer 
characteristics  (as  discussed  in  this  chapter). 

A.  EFFECT  OF  GRAVITY 

As  previously  noted  (Section  I. A),  some  improvement  in  condensation  heat  transfer 
is  expected  for  an  elliptical  tube  in  quiescent  vapor  compared  to  a  circular  one  since 
more  of  the  surface  is  aligned  with  gravity.  By  increasing  the  "effect"  of  gravity  on  the 
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condensate  film,  the  mean  condensate  velocity  at  a  given  stream  wise  location  is 
increased  which  results  in  a  thinning  of  the  condensate  film  and  increased  heat 
transfer.  The  gravity  function,  f  1(4>),  defined  in  Equation  (11-21)  is  shown  in  Figure  IV- 
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As  the  elliptical  tube  approaches  a  flat  plate  configuration  ( k  -0),  the  effect  of  gravity 
increases  to  sweep  the  condensate  film  along  the  tube  surface.  The  effect  of  gravity  on 
the  local  film  thickness  is  shown  in  Figure  IV-2for  two  eccentricities.  As  can  be  seen, 
the  condensate  film  is  thinned  over  the  front  and  rear  portion  of  the  tube  compared  to 
a  circular  tube  but  is  slightly  thicker  in  the  middle  region.  The  thickness  of  the  film 
is  controlled  by  the  rate  of  condensation  (which  is  dependent  on  condensate  film 
thickness)  and  the  condensate  velocity  (which  is  also  dependent  on  condensate  film 
thickness  and  gravity).  Near  the  top  of  the  elliptical  tube,  the  larger  gravity  component 
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increases  the  condensate  velocity  relative  to  that  of  a  circular  tube  resulting  in  a 
thinner  condensate  film.  This  thinner  film,  however,  results  in  an  increased 
condensation  rate  which  tends  to  thicken  the  film  further  downstream.  Since  the 


elliptical  tube  has  a  greater  condensation  rate  near  the  top  of  the  tube,  the  film  thickens 


Figure  FV-2.  Local  Film  Thickness  for  Varying  Eccentricity  ( k ). 


more  rapidly  compared  to  that  on  a  circular  tube  causing  the  relatively  thicker  film  in 
the  middle  region  as  seen  in  the  figure.  At  90°,  the  gravity  component  is  the  same  as  on 
a  circular  tube  and  so  has  no  relative  thinning  effect.  Over  the  rear  half  of  the  tube, 
the  reduced  condensation  rate  (as  a  result  of  this  thicker  condensate  film)  and  increased 
velocity  due  to  the  larger  gravity  component  relative  to  a  circular  tube  results  in 
continued  thickening  of  the  film  on  the  elliptical  tube,  but  at  a  slower  rate  than  that 
on  a  circular  tube.  Eventually,  the  elliptical  tube  film  thickness  is  again  thinner  than 
on  a  circular  tube  for  the  same  dimensionless  streamwise  length.  The  condensate  film 
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thickness  is  therefore  determined  by  a  balance  of  the  effects  of  two  effects:  increased 

condensate  velocity  (which  tends  to  thin  the  film)  and  increased  condensation  rate 

(which  tends  to  thicken  the  film). 

The  overall  effect  on  the  mean  heat-transfer  coefficient  (given  as  a  mean  Nusselt 
number)  is  shown  in  Figure  IV-3.  To  compare  the  elliptical  tube  with  a  circular  tube, 

the  leading  coefficient  of  the  Nusselt  solution  (Equation  I- 3b)  is  plotted  against 


eccentricity.  Solutions  for  a  circular  cylinder  (k  =  1.0)  and  a  vertical  flat  plate  (k  -  0, 
L  -  2a)  agree  well  with  the  results  of  Nusselt  [5].  For  the  horizontal  flat  plate  ( k  -  »), 
the  heat- transfer  coefficient  approaches  zero  as  would  be  predicted  by  the  Nusselt 
model.  For  the  practical  range  of  eccentricities  (0.3  <  k  <  0.6),  it  can  be  seen  that  the 
effect  of  placing  more  of  the  tube  surface  in  the  direction  of  gravity  is  to  increase  the 
mean  heat- transfer  coefficient  by  7%  for  k  =  0.6  and  13%  for  k  =  0.3.  Conversely,  by 
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placing  more  of  the  tube  surface  perpendicular  to  the  direction  of  gravity  (i.e.. 
approaching  a  horizontal  flat  plate),  the  mean  heat-transfer  coefficient  decreases  as 
expected.  These  results  are  in  agreement  with  earlier  elliptical  tube  studies  [32,33]. 


B.  EFFECT  OF  VAPOR  VELOCITY 

1.  Asymptotic  Interfacial  Shear  Stress  Approximation 

The  streamlined  shape  of  the  elliptical  tube  has  an  effect  on  the  vapor  flow 
over  the  tube.  The  vapor  velocity  function,  f2(<j>),  defined  by  Equation  (II -47)  is  shown 
in  Figure  IV-4  versus  dimensionless  streamwise  length  x  .  As  seen  in  the  figure,  the 
elliptical  tube  experiences  higher  vapor  velocities  than  a  circular  tube  at  the  front  and 
rear  of  the  tube  but  a  lower  vapor  velocity  in  the  central  region.  Using  a  balance  of  the 


factors  which  thicken  and  thin  the  condensate  film,  for  the  case  of  pure  forced 
convection  ( F  -0,  no  gravity  effect),  the  higher  relative  vapor  velocity  at  the  top  of  the 
elliptical  tube  results  in  a  larger  interfacial  shear  stress  and  hence  a  thinner  condensate 
film  and  larger  condensation  rate  compared  to  a  circular  tube.  With  more  condensate 
flowing  into  the  middle  section  of  the  elliptical  tube  and  a  lower  relative  vapor 
velocity,  the  condensate  film  thickens  more  rapidly  than  in  the  case  of  a  circular  tube. 
As  the  condensate  flows  over  the  rear  of  the  tube,  the  decreased  condensation  rate  (due 
to  the  thicker  film  in  the  middle  region)  and  increased  vapor  shear  results  in  thinner 
film  compared  to  the  circular  tube.  Figure  IV -5  shows  these  effects  in  detail  for  a 


Figure  IV-5.  Effect  of  Vapor  Velocity  on  Local  Film  Thickness  for  a  Circular 
Tube  and  Elliptical  Tube  ( k  =  0.6). 


circular  tube  ( k  =  1)  and  an  elliptical  tube  (k  =  0.6)  using  the  asymptotic  interfacial 
shear  stress  (Sbekriladze-Cinmelauri)  approximation  and  the  two-phase  boundary  layer 
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shear  stress  (Fujii  et  al.)  approximation  (to  be  discussed  later).  Whereas  the  gravity 
effect  (f](<t>))  described  in  Section  IY.A  was  always  larger  for  the  elliptical  tube  (see 
Figure  IV'- 1),  the  vapor  shear  effect  f2(4>)  is  larger  or  smaller  depending  on  the  section 
of  the  tube  to  be  analyzed  (see  Figure  IV-4).  The  overall  effect  of  vapor  shear  using  the 
Shekriladze-Gomelauri  [21]  model  is  to  reduce  (by  about  2%  in  the  practical  range  of 
eccentricities)  the  elliptical  tube  mean  heat-transfer  coefficient  relative  to  an 
equivalent  circular  tube. 

In  the  mixed  convection  region,  the  mean  heat-transfer  coefficient  of  an  elliptical 
tube  is  increased  or  decreased  compared  to  a  circular  tube  depending  on  the  relative 
magnitudes  of  the  vapor  shear  and  gravity  effects.  This  relative  magnitude  is  measured 
by  dimensionless  parameter,  F.  For  large  F,  gravity  is  dominant  and  the  heat- transfer 
coefficient  trend  is  described  in  Section  IV. A  .  For  small  F,  vapor  shear  is  dominant 
and  the  trends  described  in  this  section  determine  the  heat  transfer.  Figure  IV-6sbows 


Figure  IV-6.  Effects  of  Vapor  Velocity  on  Mean  Nu  Using  the  Asymptotic 
Interfacial  Shear  Stress  Approximation  Analysis  Method. 
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the  effect  of  F  on  the  mean  heat-transfer  coefficient  for  varying  eccentricity.  As  F  - 
<*>,  corresponding  to  free  convection  condensation,  solutions  for  elliptical  and  circular 
tubes  match  those  reported  in  Section  IV.A.  For  the  case  of  a  circular  tube  for  all 
values  of  F,  the  present  numerical  solutions  are  within  0.4%  of  an  equation  by  Rose  [28]: 


Nute1*  -  -0-9  +  -P-728  F1/2 

(  1  +  3.44  Fw  +  F  )1/4 


(IV- 2) 


2.  Two-Phase  Boundary- Layer  Shear  Stress  Approximation 

Vapor  boundary -layer  separation  is  not  predicted  by  the  asymptotic  shear 
stress  approximation  since  the  interfacial  shear  is  based  on  potential  flow  outside  the 
vapor  boundary-layer  which  is  always  positive.  Solution  of  the  two-phase  boundary- 
layer  equations  for  the  condensate  and  vapor,  with  matched  shear  stress  at  the  interface, 
allows  determination  of  the  vapor  boundary -layer  separation  point  and  its  effect  on  the 
condensation  heat  transfer.  The  analysis  also  introduces  a  dimensionless  parameter  G 
which  is  proportional  to  the  condensation  rate.  Figure  IV  -5  compares  the  asymptotic 
shear  stress  approximation  with  the  two -phase  boundary  layer  shear  stress  approxima¬ 
tion  for  a  large  condensation  rate  (large  G)  for  a  circular  and  elliptical  tube  ( k  =  0.6). 
Over  the  forward  and  middle  sections  of  the  tube,  the  general  trends  are  as  described 
in  Section  IV.B.l  for  the  asymptotic  shear  stress  approximation.  Over  the  rear  of  the 
tube,  however,  as  the  separation  point  is  approached,  the  condensate  film  thickens  more 
rapidly  due  to  the  reduced  vapor  shear  effect  (r  -  0).  It  can  be  seen  that  for  the 
elliptical  tube,  the  separation  point  occurs  further  downstream  than  for  the  circular 
tube  and  thus  this  rapid  thickening  is  delayed.  By  delaying  separation,  more  of  the  tube 
surface  has  a  thinner  film.  This  serves  to  increase  the  heat-transfer  coefficient  for  an 
elliptical  tube  compared  to  a  circular  tube  if  it  were  the  only  effect  considered. 
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In  total,  however,  the  mean  heat- transfer  coefficient  is  affected  by  a 
combination  of  gravity,  vapor  shear  and  boundary-layer  separation.  Figure  IV-7 shows 
the  effect  of  F,  G  and  k  on  the  vapor  boundary- layer  separation  point  (xs  ).  As 
previously  explained  (Chapter  III),  solutions  could  not  be  obtained  for  all  possible  Fand 
G  combinations  (particularly  large  F  and  G).  In  general,  as  F  and  G  increase  (less  vapor 
shear  and  higher  condensation  rate),  or  k  decreases  (more  elliptical),  the  vapor 
separation  point  shifts  toward  the  rear  of  the  tube.  For  low  condensation  rates  (small 
G)  corresponding  to  low  vapor  suction,  vapor  separation  occurs  near  the  positions 
obtained  for  single  phase  separation  without  suction.  Little  change  occurs  in  the 
separation  point  location  as  vapor  velocity  is  varied.  For  large  condensation  rates  (large 
G),  the  vapor  boundary -layer  experiences  high  suction  which  shifts  the  separation  point 
downstream  (i.e.,  as  G  increases,  ;rs  increases).  As  vapor  velocity  is  decreased 
(increasing  F),  the  separation  point  shifts  further  downstream.  As  k  decreases  (more 
elliptical),  the  vapor  boundary- layer  separation  point  moves  further  downstream  when 
compared  to  a  circular  tube  under  all  conditions  of  F  and  G.  At  high  F  and  G,  the  vapor 
boundary- layer  separation  points  come  together  at  the  rear  of  the  tube  for  both  the 
elliptical  and  circular  tube. 

Figure  IV-8shows  the  influence  that  the  above  effects  (vapor  shear,  gravity, 
condensation  rate  and  eccentricity)  have  on  the  mean  heat-transfer  coefficient.  The 
solutions  for  k  =  1  agree  with  the  results  of  Fujii  et  al.  [22]  as  corrected  by  Lee  and  Rose 
[25].  Similar  trends  were  obtained  for  the  elliptical  tubes.  For  large  values  of 
condensation  rate  (large  G),  the  two- phase  boundary -layer  shear  stress  analysis  agrees 
closely  with  the  asymptotic  shear  stress  analysis  method  as  expected  since  the  latter 
assumes  an  infinite  condensation  rate.  For  decreasing  k,  the  results  for  an  elliptical 
tube  tend  to  show  a  reduction  in  heat  transfer  at  low  F  (compared  to  a  circular  tube) 


66 


Figure  IV- 7.  Vapor  Separation  Point  Location  with  Varying  F,  G  and  k. 
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Figure  IV-8.  Effects  of  Vapor  Velocity  and  Vapor  Separation  on  Mean  Nu  using 
the  Two-Phase  Boundary- Layer  Approximation  Analysis  Method. 


and  an  increase  in  heat  transfer  at  high  F,  as  previously  mentioned.  The  cross-over 
point  from  reduction  to  improvement  occurs  at  increasing  values  of  F  as  G  decreases. 
For  large  values  of  F,  the  results  converge  to  the  Nusselt  analysis.  This  shift  in  the 
cross-over  point  is  attributed  to  the  change  in  the  vapor  boundary- layer  separation  point 
with  varying  condensation  and  vapor  velocity  (varying  G  and  F).  As  seen  in  Figure  IV- 
7,  for  large  condensation  rate  ( G  =  5.5)  and  low  vapor  velocity  (large  F),  the  large  vapor 
suction  causes  a  delay  in  vapor  boundary- layer  separation  such  that  it  occurs  at 
essentially  the  same  streamwise  location  for  all  eccentricities.  As  velocity  increases  (F 
becomes  smaller)  the  separation  point  shifts  forward  at  different  rates  dependent  on 
eccentricity  which  is  evidenced  by  the  diverging  curves  in  Figure  IV-7c.  Since  early 
vapor  boundary-layer  separation  results  in  a  reduction  in  heat  transfer  due  to  rapid 
thickening  of  the  condensate  film,  the  circular  tube  ( k  =  l)  experiences  more  of  a 
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reduction  in  heat  transfer  compared  to  an  elliptical  tube.  The  effect  of  delaying 
separation  is  to  cause  a  thinning  in  the  mean  film  thickness  for  the  tube.  In  the 
balance,  where  the  streamwise  component  of  gravity  and  the  shift  in  vapor  boundary- 
layer  separation  are  positive  effects  (i.e.,  both  contribute  to  a  relative  thinning  of  the 
elliptical  tube  condensate  film  compared  to  a  circular  tube)  and  vapor  shear  is  a 
negative  effect  (i.e.,  lower  shear  contributes  to  a  relative  thickening  of  the  elliptical 
tube  condensate  film  compared  to  a  circular  tube),  delays  in  vapor  separation  move  the 
cross-over  point  from  gravity  dominant  to  shear  dominant  heat  transfer.  For  the  case 
of  low  G  in  Figure  IV-7a,  there  is  no  noticeable  shifting  of  the  vapor  boundary- layer 
separation  point  as  F  is  varied.  Thus  the  vapor  boundary-layer  separation  point 
contribution  does  not  change  with  vapor  velocity  (F)  to  the  balance  and  the  cross-over 
point  occurs  at  lower  velocities  (larger  F).  The  overall  effect  of  vapor  shear  using  this 
analysis  is  typically  a  slight  reduction  (<  2%)  in  the  heat-transfer  coefficient  for  an 
elliptical  tube  compared  to  an  equivalent  surface  area  circular  tube,  dependent  on  the 
magnitude  of  the  condensation  rate  parameter,  G. 

Figures  IV-9and  IV-lOshow  the  effects  of  vapor  boundary-layer  separation 
on  the  local  film  thickness  and  heat- transfer  coefficient  for  a  given  dimensionless 
streamwise  distance.  For  low  F  and  G,  the  separation  point  occurs  relatively  early,  but 
is  significantly  delayed  when  using  an  elliptical  tube.  As  F  increases  (vapor  velocity 
decreases),  the  separation  points  remain  at  their  same  respective  streamwise  locations. 
The  areas  under  the  curves  in  Figure  IV- 10  represent  the  total  heat  transfer.  Bearing 
this  in  mind,  the  comparison  between  elliptical  and  circular  tubes  is  made  a  little  easier. 
At  high  F  (for  given  G),  the  area  under  the  curve  for  an  elliptical  tube  exceeds  that  for 
a  circular  tube  while  at  low  F,  the  opposite  is  true.  However,  for  a  given  eccentricity, 
the  values  of  jcs  are  the  same  for  varying  F  (as  noted  earlier). 
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As  G  increases  (increasing  condensation  rate  and  suction),  vapor  boundary-layer 
separation  is  further  delayed  (compared  to  low  G)  and  is  virtually  eliminated  at  high 
F.  Separation  is  always  delayed  for  an  elliptical  tube  compared  to  a  circular  tube. 
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Figure  IV-9cand  d.  Dimensionless  Local  Film  Thickness  tor  Varying  G,  F  and  k 
from  the  Two-Phase  Boundary  Layer  Shear  Stress  Method  of  Analysis. 
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Figure  10  a  and  b.  Local  Nu  for  Varying  G,  F  and  k  from  the  Two-Phase  Boundary 
Layer  Shear  Stress  Analysis  Method. 
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Figure  10  c  and  d.  Local  Nu  for  Varying  G.  F  and  k  from  the  Two-Phase  Boundary 
Layer  Shear  Stress  Analysis  Method. 
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C.  EFFECT  OF  SURFACE  TENSION 


As  noted  in  Section  II. E,  a  condition  exists  for  which  d6  /d<f>  becomes  infinite  at 
some  critical  angle  <j>c.  What  is  actually  happening  can  be  discerned  from  Figure  IV- II 
which  shows  the  surface  tension  function  f4(<|>)  versus  .t  for  varying  eccentricity. 
Surface  tension  causes  a  favorable  pressure  gradient  over  the  front  half  of  the  elliptical 
tube  and  an  adverse  pressure  gradient  over  the  back  half.  The  severity  of  the  pressure 
gradient  is  localized  to  small  regions  at  the  top  and  bottom  of  the  elliptical  tube  where 


Figure  IV -11.  Surface  Tension  Function,  f4(<J>)  for  Varying  Eccentricity. 


change  in  surface  curvature  is  most  severe.  Therefore,  it  is  most  significant  for  small 
values  of  eccentricity.  The  above  analysis  explains  the  reason  for  condensate  flow 
instability  discussed  in  Section  II. E  and  predicted  by  Equation  (11-90). 
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Further  discussion  of  surface  tension  is  restricted  to  those  cases  in  which  a 
solution  could  be  obtained  over  at  least  90%  of  the  tube  surface.  This  restriction  was 
arbitrarily  chosen  based  on  the  small  amount  of  heat  transfer  which  occurs  on  the  lower 
part  of  the  tube  and  would  therefore  have  minimal  impact  on  the  accuracy  of  the 
results.  Figure  IV- 12  shows  the  effects  of  surface  tension  on  the  mean  heat-transfer 
coefficient  for  k  =  0.6  and  0.4  .  The  Bond  number,  Bo,  gives  the  relative  effect  of 
inertia  to  surface  tension.  Values  of  1  IBo  of  0.01  (typical  for  steam)  and  0.001  (typical 
for  refrigerants)  as  well  as  1  /Bo  =  0  have  been  shown  on  the  figure.  A  number  of 
interesting  points  can  be  highlighted.  Firstly,  the  surface  tension  effect  is  much  smaller 
for  highly  wetting  refrigerants,  as  expected.  Secondly,  inclusion  of  surface  tension 
leads  to  a  small  decrease  in  the  mean  heat-transfer  coefficient  (<  2%)  over  the  whole 
range  of  F  for  the  practical  range  of  eccentricities,  suggesting  that  any  thinning  of  the 
condensate  film  over  the  top  half  of  the  tube  is  more  than  offset  by  a  thickening  over 
the  lower  half.  As  eccentricity  decreases,  this  discrepancy  is  accentuated.  Finally, 
surface  tension  effects  are  felt  more  in  the  free  convection  region  (high  F)  than  in  the 
forced  convection  region  (low  F)  as  film  thickness  becomes  dominated  by  vapor  shear. 

Krupiczka  [29]  used  the  curvature  of  the  condensate  film  surface  to  analyze 
surface  tension  on  a  horizontal  circular  cylinder.  He  determined  that  the  effect  of 
surface  tension  was  only  significant  over  the  bottom  portion  of  the  tube.  This  fact 
implies  that  since  the  film  is  very  thin  over  the  top  of  the  tube,  its  curvature  is  very 
close  to  that  of  the  tube  surface  itself,  which,  for  a  circular  tube  results  in  no  surface 
tension  effect  (constant  radius  of  curvature).  As  the  condensate  flows  around  the  lower 
half  of  the  circular  tube,  the  film  curvature  no  longer  follows  the  tube  surface.  Rather, 
th^  radius  of  curvature  is  increasing  resulting  in  a  small  improvement  in  the  heat- 
transfer  coefficient.  A  similar  analysis  can  be  applied  to  the  elliptical  tube.  Over  the 
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top  half  of  the  tube  where  the  condensate  film  is  thinnest,  modeling  the  surface  tension 
effect  using  the  curvature  of  the  tube  surface  is  valid.  However,  the  adverse  pressure 
gradient  over  the  lower  half  of  the  tube  is  probably  not  as  severe  as  seen  in  Figure  IV- 
11.  Over  this  portion  of  the  tube,  the  film  is  thicker  and  its  decreasing  surface 
curvature  is  not  as  severe  as  the  tube  surface  curvature.  The  conclusion  drawn  from 
this  analysis  is  that  surface  tension  does  not  have  a  significant  effect  on  the  mean  heat- 
transfer  coefficient  for  an  elliptical  tube. 

D.  EFFECT  OF  PRESSURE  GRADIENT 

The  effect  of  pressure  gradient  was  analyzed  using  the  asymptotic  shear  stress 
approximation  (similar  to  Rose  [28]  for  a  circular  tube).  In  this  analysis,  the  pressure 
gradient  in  the  condensate  film  was  assumed  to  be  due  to  the  pressure  gradient 
impressed  on  the  condensate  by  vapor  potential  flow.  As  in  the  case  of  surface  tension, 
there  is  a  relationship  between  parameters  P  and  F  (Equation  II -89)  for  which  a  solution 
could  not  be  obtained  for  the  entire  tube.  Figure  IV- 13  shows  the  pressure  gradient 
function  versus  jc*  for  varying  k  which  sets  up  a  favorable  pressure  gradient  over  the 
top  of  the  tube  and  an  adverse  pressure  gradient  over  the  lower  half  of  the  tube.  The 
effect  of  eccentricity  is  to  shift  the  point  of  maximum  pressure  gradient  to  the  front 
and  rear  of  the  tube.  Thus,  <}>c  shifts  downstream  for  an  elliptical  tube  compared  to  a 
circular  tube  for  similar  P  and  F.  This  is  a  result  of  the  more  streamlined  shape  and  the 
smaller  resultant  pressure  drop  over  the  streamwise  length  of  an  elliptical  tube. 

As  in  the  analysis  of  surface  tension,  discussion  is  limited  to  those  cases  for  which 
solutions  could  be  obtained  over  the  entire  tube.  When  pressure  gradient  was  included 
in  the  momentum  equation  of  the  condensate  film,  solutions  agreed  with  the  results  of 
Rose  [28]  for  a  circular  tube,  yielding  a  maximum  5%  decrease  in  the  heat- transfer 
coefficient  over  the  entire  range  of  F  compared  to  the  case  with  pressure  gradient 
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Figure  IV- 13.  Pressure  Gradient  Function,  f3(<J>),  for  Varying  Eccentricity. 


neglected.  Figure  IV-14  shows  the  mean  heat-transfer  coefficient  for  a  circular  tube 
(k  =  1)  for  varying  F  and  P.  The  favorable  pressure  gradient  increases  the  mean 
condensate  velocity  resulting  in  a  thinner  film  and  an  increased  mean  heal- transfer 
coefficient  over  the  lop  of  the  tube.  This  increase  is  offset  by  the  adverse  pressure 
gradient  which  slows  the  condensate  film  velocity  over  the  lower  half  of  the  tube, 
thickens  the  film  and  decreases  the  mean  heat- transfer  coefficient.  For  an  eccentricity 
of  0.6,  Figure  IV-15  compares  the  effect  of  pressure  gradient  for  an  elliptical  tube  with 
those  for  a  circular  tube  for  varying  F  and  P.  For  both  the  circular  and  elliptical  tubes, 
the  effect  of  pressure  gradient  provides  a  slight  reduction  in  condensation  heat  transfer. 
Figure  IV-15  indicates  that  the  reduction  in  heat  transfer  was  less  for  the  elliptical  tube 
compared  to  the  circular  tube  as  a  result  of  the  more  favorable  potential  velocity 
distribution  around  the  ellipse. 
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E. 


EFFECT  OF  VAPOR  PRESSURE  DROP 


The  previous  results  indicate  that  vapor  shear  decreases  the  heat-transfer 
coefficient  for  an  elliptical  tube  when  compared  to  a  circular  tube.  However,  a 
consideration  not  yet  discussed  is  that  the  vapor  pressure  drop  across  an  elliptical  tube 
is  significantly  smaller  than  a  circular  tube  for  the  same  free-stream  velocity. 
Alternatively,  for  a  given  pressure  drop,  the  free-stream  velocity  is  greater  for  the 
elliptical  tube.  Since  very  little  data  is  available  for  two- phase  drag  coefficients  around 
tubes,  single- phase  data  was  used  where  the  drag  coefficient  for  the  tube  is  due  entirely 
to  form  drag.  Values  of  CD  for  a  circular  ( k  =  1)  and  an  elliptical  tube  ( k  =  0.5)  were 
estimated  to  be  approximately  1.2  and  0.6  respectively  (White  [44]).  For  steam  with  Tsat= 
60°C  condensing  on  a  horizontal  tube  with  Twa„=  40°C  and  U„=  25  m/s  for  the  circular 
tube,  the  corresponding  vapor  velocity  for  the  elliptical  tube,  for  the  same  pressure 
drop,  was  calculated  to  be  U„=  35  m/s.  The  resulting  values  of  F  are  therefore  0.0257 
and  0.0132  for  the  circular  and  elliptical  tubes  respectively.  Using  the  asymptotic  shear 
stress  assumption  for  both  values  of  F  gives  an  increase  in  the  heat  transfer  for  an 
elliptical  tube  of  16.3%  when  compared  to  a  circular  tube.  Using  the  two- phase 
boundary -layer  shear  stress  assumption  gives  a  corresponding  increase  of  17.1%. 
Therefore,  though  the  effect  of  vapor  shear  alone  results  in  a  reduction  of  the  heat- 
transfer  coefficient  (as  discussed  in  Section  1V.B),  when  the  streamlined  shape  of  the 
elliptical  tube  and  its  effect  on  the  vapor  pressure  drop  is  taken  into  account  (allowing 
for  higher  U„  for  a  given  pressure  drop),  there  is  an  increase  in  the  heat- transfer 
coefficient  when  compared  to  a  circular  cylinder. 
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F.  INSIDE  HEAT  TRANSFER  COEFFICIENT 


So  far  only  the  outside  heat- transfer  coefficient  has  been  discussed.  In  real 
condensers,  it  is  the  overall  heat- transfer  coefficient  which  is  important  and  so 
consideration  must  also  be  given  to  the  single- phase  convective  heat  transfer  occurring 
inside  the  elliptical  tube.  A  preliminary  survey  of  the  literature  (Incropera  and  Dewitt 
[45])  shows  that  for  turbulent  flow  conditions  where  the  Prandtl  number  is  greater  than 
0.5  (which  is  true  for  most  practical  condensers),  correlations  developed  for  circular 
tubes  may  be  used  with  good  accuracy  to  approximate  the  inside  heat-transfer 
coefficient  for  an  elliptical  tube  if  the  hydraulic  diameter  is  used  in  place  of  the 
circular  diameter.  The  hydraulic  diameter  is  given  by 

(IV-3) 


4  Ar 


P 


where  Ac  is  the  cross-sectional  area  of  the  tube  and  p  is  the  inside  perimeter  of  the  tube. 

The  perimeter  for  an  elliptical  tube  and  an  equivalent  surface  area  circular  tube  is  the 

same.  The  cross-sectional  area  of  an  elliptical  tube  is  given  by 

Ac  =  nab  (IV -4) 

and  is  smaller  than  an  equivalent  surface  area  circular  tube.  Inside  heat  transfer 
correlations  for  turbulent  flow  heat  transfer  inside  a  circular  tube  are  typically  of  the 
form  [45] 


Nu 


' 

«  Reot 


P 


(IV -5) 


Thus,  the  inside  heat-transfer  coefficient  is  related  to  the  hydraulic  diameter  by 


,-0.2 


(IV-6) 
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As  an  example,  consider  an  elliptical  tube  with  k  =  0.5  and  a  major  axis  with  a  =  0.01 
m.  The  equivalent  circular  tube  would  have  a  diameter  of  0.0154  m.  The  ratio  of  the 
heat  transfer  coefficient  of  the  elliptical  tube  to  the  circular  tube  is  given  by 

^  elliptical 
circle 


elliptical 

-0.2 

^ c\uipncai  | 

^/)circU 

~ 

circle 

(I  V  -  7) 


thus, 


a elliptical 
“circle 


ab 


(DJ2)2 


-02 


=  1.035 


(IV- 8) 


For  this  simple  analysis,  it  can  be  seen  that  the  inside  heat-transfer  coefficient  of  the 
elliptical  tube  is  3.5%  greater  than  the  equivalent  surface  area  circular  tube.  This 
result,  when  added  to  the  increase  in  the  outside  heat-transfer  coefficient  (for  a  given 
pressure  drop),  would  increase  the  overall  heat-transfer  coefficient  of  the  condenser. 
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V.  CONCLUSIONS 


Analyses  of  laminar  film  condensation  on  a  horizontal  elliptical  tube  were 
conducted  under  conditions  of  free  and  forced  convection,  the  latter  using  both  an 
asymptotic  interfacial  shear  stress  approximation  as  well  as  solution  of  the  two- phase 
boundary -layer  equations.  For  the  asymptotic  shear  stress  approximation,  the  effects 
of  surface  tension  and  pressure  gradient  were  also  included.  Where  possible,  these 
analyses  have  been  validated  against  existing  solutions  for  laminar  film  condensation 
on  a  horizontal  circular  tube  and  a  vertical  flat  plate. 

Whether  or  not  condensation  heat  transfer  is  improved  for  an  elliptical  tube 
compared  to  an  equivalent  surface  area  circular  tube  is  dependent  on  a  balance  of  the 
effects  of  gravity,  vapor  shear  and  vapor  boundary -layer  separation  and  the  influence 
these  have  on  the  thickness  of  the  condensate  film.  Under  quiescent  vapor  conditions 
(no  vapor  shear  or  boundary-layer  separation),  gravity  causes  an  increase  (-10%)  in  the 
heat  transfer  on  an  elliptical  tube  compared  to  a  circular  tube  of  the  same  surface  area 
due  to  an  increase  in  "effective"  gravity.  Under  conditions  of  forced  convection,  both 
shear  stress  approximations  indicate  a  small  decrease  (<  2%)  in  performance  due  to  the 
reduction  in  interfacial  shear  as  a  result  of  the  better  streamlined  shape  of  the  elliptical 
tube  when  compared  to  the  circular  tube  in  the  same  free  stream  vapor  flow.  However, 
when  pressure  drop  effects  are  also  considered,  the  higher  allowable  vapor  velocity  over 
an  elliptical  tube  (for  the  same  pressure  drop  as  a  circular  tube)  results  in  an  increase 
in  the  mean  heat-transfer  coefficient  of  15-20%.  For  conditions  of  mixed  convection, 
the  condensate  film  thickness  is  controlled  by  both  gravity  and  vapor  shear  effects.  The 
cross-over  between  heat  transfer  improvement  (gravity  dominant  flow)  and  heat 


transfer  reduction  (vapor  shear  dominant  flow)  for  the  elliptical  tube  is  dependent  on 
the  vapor  boundary- layer  separation  point.  For  large  condensation  rates  (high  vapor 
suction),  vapor  boundary-layer  separation  is  delayed  and  the  cross-over  point  shifts  to 
higher  velocities  (smaller  F).  For  low  condensation  rates  (low  vapor  suction),  movement 
of  the  vapor  boundary -layer  separation  point  is  minimal,  resulting  in  a  minor  effect  on 
heat  transfer  and  a  cross-over  point  that  shifts  toward  lower  velocities  (large  F). 
Pressure  gradient  and  surface  tension  each  lead  to  a  small  decrease  in  the  mean  heat- 
transfer  coefficient  (<  2%)  for  an  elliptical  tube. 

In  general,  therefore,  the  outside  heat-transfer  coefficient  of  a  condenser  tube  is 
improved  by  using  an  elliptical  tube  geometry.  Approximation  of  the  single -phase 
inside  heat-transfer  coefficient  using  a  hydraulic  diameter  and  a  Seider-Tate  type 
correlation  indicates  that  the  elliptical  tube  also  has  better  inside  heat  transfer 
performance  compared  to  an  equivalent  surface  area  circular  tube.  Thus,  the  overall 
heat- transfer  coefficient  should  be  improved. 
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VI.  RECOMMENDATIONS 


This  study  represents  an  initial  step  in  determining  if  a  shell-and-tube  condenser 
using  horizontal  elliptical  tubes  has  better  heat  transfer  than  a  comparable  circular  tube 
bundle.  The  effects  of  multi-tube  arrangements  (i.e.,  the  influence  of  adjacent  tubes 
and  the  spacing  between  tubes)  on  vapor  flow  characteristics  needs  to  be  evaluated  to 
determine  if  further  improvements  in  the  overall  heat-transfer  coefficient  can  be 
achieved.  It  is  suspected  that  the  interfacial  shear  of  the  condensate  film  on  elliptical 
tubes  will  be  further  reduced  but  the  corresponding  reduction  in  vapor  pressure  drop 
across  the  condenser  will  result  in  a  net  gain  in  heat  transfer  compared  to  a  circular 
tube  bundle.  An  analysis  similar  to  that  conducted  by  Aoune  and  Burnside  [46,  47]  can 
be  used  to  evaluate  multi -tube  arrangements. 

Further  studies  should  be  conducted  on  a  single  elliptical  tube  to  better  determine 
the  effect  of  surface  tension  on  heat  transfer.  The  surface  tension  model  used  in  this 
study  was  based  only  on  the  curvature  of  the  tube  surface  and  does  not  take  into 
account  the  relatively  thicker  film  over  the  lower  half  of  the  tube.  The  curvature  of 
the  actual  condensate  film  surface  should  be  analyzed  in  a  manner  similar  to  Krupiczka 
[29].  This  ana. » sis,  together  with  the  asymptotic  shear  stress  assumption,  would  involve 
the  solution  of  a  second  order  ODE. 

Analysis  of  the  effect  of  offsetting  the  elliptical  tube  at  an  angle,  a,  with  respect 
to  vapor  flow  and  gravity  must  be  conducted  to  determine  its  effect  on  the  heat- 
transfer  performance  of  the  tube.  This  is  important  since  in  real  condensers,  vapor  will 
approach  the  tube  at  various  angles.  Additionally,  since  the  elliptical  tube  is  a 
streamlined  body,  vapor  flow  in  a  direction  other  than  the  direction  of  the  major  axis 
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may  result  in  a  lifting  force  on  the  tube  which  may  result  in  tube  vibration  and  wear. 
These  effects  would  have  to  be  studied  and  methods  to  minimize  any  resultant  vibration 
must  found.  The  problems  associated  with  joining  elliptical  tubes  to  the  condenser  tube 
sheet  and  baffle  plates  also  needs  to  be  addressed. 

Finally,  experiments  must  be  conducted  on  elliptical  lubes  to  validate  the 
predictions  of  these  theoretical  models. 
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APPENDIX  A  -  DERIVATION  OF  TANGENTIAL  VAPOR  VELOCITY  OVER  AN 


ELLIPTICAL  TUBE  FROM  POTENTIAL  FLOW  THEORY 


The  basic  technique  used  to  determine  fluid  flow  velocity  over  a  body  is  to 
employ  a  conformal  transformation  to  translate  the  body  from  one  complex  plane  into 
another  complex  plane.  From  this,  the  velocity  of  the  transformed  body  can  be  readily 
computed  using  potential  theory.  For  this  particular  application,  the  fluid  velocity 
along  the  surface  of  a  circular  cylinder  is  determined  and  then  transformed  into  the 
complex  plane  in  which  the  cylinder  has  an  elliptical  shape. 

Consider  a  circular  cylinder  in  the  complex  plan*,  Zj,  as  shown  in  Figure  A-l.The 
surface  of  the  circle  is  defined  by 


z,  -  c  e16  (a) 


(A  - 1) 


-  y]  -  r*  (b) 

The  conformal  transform  to  translate  the  circular  body  from  the  zx  plane  to  an  elliptical 
body  in  the  z2  plane  is  given  by 


h  =  *1 


*1 


(A-2) 


thus 


*2  =  <*1  +  =  *,  1  -  \ 
xi  +  Vi  c1) 


Vi 


(A-3) 


The  real  and  imaginary  coordinates  in  the  z2  plane  are  given  by 


x2  =  *, 


and  y2  =  y, 


1  +  — 
c2 


(A -4) 
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where  Uj  and  Vj  are  the  fluid  velocity  components  in  cartesian  coordinates.  The  fluid 


velocity  in  the  z,  plane  (for  the  ellipse)  is  found  by  the  chain  rule,  thus, 

dw  dwdz  i 

—  = - =  u,  -  IV, 

dzxdz1 


(A -8) 


Combining  Equations  (A-la),(A-2)and  (A-7)  results  in 

dw  U.  c2  (e‘ &  - -  e11"2) 
c2  ei26  +  1 


(A -9) 


The  magnitude  of  the  fluid  velocity  is  found  by  multiplying  the  complex  velocity  by 


its  conjugate.  Using  Euler’s  formula, 

e‘e  =  cos 0  +  i  sin© 


(A- 10) 


and  with  some  algebraic  manipulation,  the  potential  velocity  of  the  fluid  at  the  surface 


of  the  ellipse  is  given  by 


U2  =  =  U~  c*  (2  *  2cos20) 

8  dzj  rfej  c*  +  l  +  2c2  cos26 


(A- 11) 


where  the  bar  denotes  the  complex  conjugate.  Equation  (A-ll)is  further  simplified 


using  the  trigonometric  identities 


cos26  =  cos2©  -  sin2© 
cos23  =  2cos2©  -  1 


(A- 12) 


and  the  relationship  between  c,  a  and  b. 


4 c2  =  (a+bj2 


(A- 13) 


Thus  the  potential  velocity  is  given  by 


U  =  U  U+*)  cosQ 


/cos2©  +  ifc2sin26 


(A- 14) 


The  angle,  0,is  measured  from  the  xx  (horizontal)  axis.  The  relationship  between  ©and 
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4>,  the  angle  measured  from  the  _V[  (vertical)  axis,  is  given  by 

0  =  -  -  4,  (A- 15) 

2 

Note  that  <f»  is  an  angle  defined  in  the  zl  plane  for  the  circle  and  is  the  parametric  angle 
defined  in  Section  II.A  .  The  potential  velocity  is 

=  um  (l+k)  sin4>  =  Um 
y/sin24>  +  t2cos24> 


thus  given  by 
(1+*) 

^1  +  k2  cot2^ 


(A- 16) 
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APPENDIX  B  -  COMPUTER  CODE  FOR  ASYMPTOTIC  SHEAR  STRESS  MODEL 


*  PROGRAM  MIXED1 

SXKXX***XKXXXXXXXXKXKXXKKKXXXXXKXXXXXKXKXXXKXXXXKX*XXXXKKXXXXXXX*XXKX* 

X 

*  THIS  PROGRAM  NUMERICALLY  SOLVES  THE  PROBLEM  OF  MIXED 

X  CONVECTION  CONDENSATION  WITH  SURFACE  TENSION  ON  AN  ELLIPTICAL 
x  CYLINDER.  THE  PROBLEM  FORMULATION  USES  THE  ASYMPTOTIC  SHEAR 
x  STRESS  APPROXIMATION  ANALYSIS  METHOD. 

X 

x  THE  PROGRAM  USES  AN  AUTHOR  DEVELOPED  FORWARD  STEPPING 
x  ALGORITHM  BASED  ON  A  4TH  ORDER  ADAMS  PREDICTOR-CORRECTOR 

X  METHOD  TO  SOLVE  A  1ST  ORDER  ORDINARY  DIFFERENTIAL  EQUATION, 

x 

X  THIS  PROGRAM  CAN  BE  RUN  WITH  DIMENSIONLESS  PARAMETERS 
X  (F,G,P,  BO)  SPECIFIED  OR  WHERE  THE  PARAMETERS  ARE  CALCULATED 
X  BASED  ON  SPECIFIED  FLUID  PROPERTIES.  MODIFICATIONS  TO  THE 
x  PROGRAM  (REMOVAL  OF  UNNECCESSARY  STEPS  BY  "COMMENTING"  THE 

x  STATEMENT  OUT  OF  THE  PROGRAM)  IS  REQUIRED.  SPECIFIED  FLUID 

x  PROPERTIES  INCLUDE  T  SAT,  T  WALL,  AND  UINF. 

X 

x  NUMERICAL  RESULTS  ARE  DIRECTED  TO  AN  EXTERNAL  DATA  FILE, 
x 

X  MAJOR  VARIABLES  USED  AREt 
x 


X 

THETA 

PARAMETRIC  ANGLE 

X 

X 

STREAMWISE  LENGTH 

X 

DEFF 

EFFECTIVE  DIAMcTER 

X 

DEL 

CONDENSATE  FILM  THICKNESS 

X 

DELND, Z 

DIMENSIONLESS  FILM  THICKNESS 

X 

F, P, BO. 

DIMENSIONLESS  PARAMETERS 

X 

jarprl 

X 

UINF 

VAPOR  FREE  STREAM  VELOCITY 

X 

TS 

VAPOR  SATURATION  TEMPERATURE 

X 

TW 

WALL  SURFACE  TEMPERATURE 

X 

DENSL 

CONDENSATE  DENSITY 

X 

DENSV 

VAPOR  DENSITY 

X 

CONDL 

CONDENSATE  THERMAL  CONDUCTIVITY 

X 

VISL 

CONDENSATE  VISCOSITY 

X 

HFG 

LATENT  HEAT  OF  VAPORIZATION 

X 

NU 

NUSSELT  NUMBER 

X 

NUAVG 

MEAN  NU  NUMBER 

X 

ALPHA 

HEAT  TRANSFER  COEFFICIENT 

X 

QFLUX 

HEAT  FLUX 

X 

A 

SEMI-MAJOR  AXIS 

X 

r 

ECCENTRICITY 

X 

XFN 

PARAMETRIC  RADIUS 

X 

G 

GRAVITY  FUNCTION 

X 

FI 

VELOCITY  FUNCTION 

X 

Fll 

DERIVATIVE,  VELOCITY  FUNCTION 

X 

F2 

PRESSURE  GRADIENT  FUNCTION 

X 

F3 

SURFACE  PENSION  FUNCTION 

X 

XXKXXXXXXXXKXXXKXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

PARAMETER(PI*3. 141592654, GRAV*9 .807 ) 

REAL  Z( 0 : 180000 ), THETA( 0 j 180000 ),X(0: 180000 ),H,DEFF, XI, X2,X3, 
+NUAVG, NU, NUORE, UINF,TS , TW, TL , DENSL , VISL , TG, CONDL , HFG, B2, DELND, 
♦DEL  ,  QFLUX,  A ,  K, RE, SUM, XFN , ALPHA, ALPHA 0 , DQFLUX, F, PHI , 
+A1,A2,A3.A4,A5,A6,A7,A8,A9,A10,A11,TF,PSAT,TB,TC,TD,R,TE, DENSV, 
♦SIGML  ,  FCN ,  P ,  Cl ,  C2 ,  C3 ,  FCO ,  FC1 ,  FC2,  FC3 , ZP1 ,  ZP2,  ZP3 , ZP4 ,  BO, 

♦C4, C5, F21 , F31 , NUCOEFF 

INTEGER  N,J,JJ,IC 

COMMON  ✓GEOM/'A ,  K ,  DEFF/PARAM/F ,  BO ,  P 

EXTERNAL  FCN . XFN , G , G1 , FI , FI 1 , F2 , F3 

OPEN (20, FILE*' /'MIXEDl  OUTPUT  A1 '  .  STATUS*  'OLD' ) 

0PEN( 50 , FILE1 '/DEL 2  DATA  A1 STATUS* ' OLD' ) 
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FILE:  MIXED1 


FORTRAN  A 


***  INPUT  FLUID  PARAMETERS  *** 

PRINT*, 'ELLIPSE  OF  THE  FORM  (X/'8)**2+(Y/A)**2=l » 

PRINT*, 'INPUT  "A"  DIMENSION  CM)?' 

READ*, A 

PRINT*, 'INPUT  T  SAT  AND  T  WALL  CIO?' 

READ*, TS , TW 

PRINT*, 'INPUT  ANGLE  STEP  SIZE  (DEGREES)?' 

READ* , H 

PRINT*, 'STEP  SIZE,  H  =  ' , H 
C  PRINT*, 'INPUT  UINF?' 

C  READ*, UINF 

PRINT*, 'INPUT  K?' 

READ*, K 

PRINT*, 'INPUT  F?' 

READ*, F 

N=180/H+l 

H=PI/180.*H 

B0=0. 

P=0. 

***  DETERMINE  FLUID  PROPERTIES  ***xxx*x 

Al=15. 49217901 
A2=-5. 6783717693 
A3=l. 4597584637 
A4=13. 877000608 
5=-80. 887673591 
A6=123. 56883468 
A7=-188. 321212064 
A8=660. 91763485 
A9=-1382. 4740091 
A10=1300. 1040184 
All=-449. 39571976 
TF=TS/1 000 . 

PSAT=1.E+6*EXP(A1+A2/TF+A3*10GCTF)+TF*A4+A5XTF**2+A6*TF*X3+A7X 
+TFX*4+A8XTF**5+A9*TF**6+A1 0*TF**7+A1 1*TF**8 ) 

TB=1500./TS 

TC=2 . 5*LOG< 1 . -EXP(-TB) ) 

TD=.0015/(1.+.0001*TS)-.000942*EXP(TB+TC)*SQRT(1./TB)-.0004882*TB 
R=461 .51 

TE=2 . XPSAT/ ( R*TS ) 

DENSV=TE/ C 1 . +SQRT( 1 . +2 . *TD*TE) ) 

TL  =  (2.*TW+TSV3. 

DENSL=1 .  0012674-TL*(  2 . 02915E-6-TLX3 . 8333E-9 )  ) 
VISL=2.414E-5*1Q**(Z47.8/(TL-140.)) 

TG=TL/27  3 . 15 

CONDL=- . 92407+TGXC2 . 8395-TG*( 1 ,8007-TG*(  .  52577-TG* . 07344) ) ) 
HFG=3468920.-TS*(5707 .4-TS*(ll . 556 2-TS* . 0133103) ) 
SIGML=(-.0003*(TL-273.15)**2-.138*(TL-273.15)+75.6)/1000. 

*X*  DETERMINE  THETAC I ) , X( I ) , DEFF  ******** 

x*X  USE  SIMPSON'S  RULE  TO  NUMERICALLY  INTEGRATE  ALONG  SURFACE  *** 

THETAC 0)=0. 

XC0)=0. 

DO  10,  1=1, N 

THETA(I)=THETA(I-1)+H 
XI =A*XFN( THETAC 1-1 ) ) 

X2=A*XFN( ( THETAC I -1 )+THETA( I ) )* . 5) 

X3=A*XFNC  THETAC I ) ) 

X  C I )=XC 1-1 )+H/6 . XCX1+4 . *X2+X3) 

10  CONTINUE 

DEFF=2.*XCN)/PI 

***  DETERMINE  DIMENSIONLESS  PARAMETERS  AS  REQUIRED  **** 

C  F=VISL*DEFF*HFG*GRAV/(UINF**2*C0NDL*(TS-TW) ) 

UINF=SQRT( VISL*HFG*GRAV/( F*CONDL*(TS-TW) ) ) 
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FILE:  MIXED1  FORTRAN  A 


C  P=DENSVXHFGXVISL/(DENSLXCONDLx(TS-TW)) 

C  30=SIGML/(DENSLXGRAVXDEFFXX2) 

RE=DENSLXDEFFXUINF/VISL 

PRINTX, 'COMPLETED  THETA, X, DEFF* 

WRITE(2Q,5000)  'SAT.  TEMP.  (K)  =  ',TS,'WALL  TEMP.  =  ',TW, 

+  'l/BO  =  ' , BO,  'P  =  ' , P , ' F  =  ' , F, 'UINF  (M/S)  =  ' , UINF, 'REL  =  *,RE 
WRITEC20 , 5010 )  '"A"  LENGTH  CM)  =  ',A,"*K"  ECCENTRICITY  = 

+ ' DEFF  (M)  =  ' , DEFF 
WRITEC20, 5020) 

xxx  DETERMINE  INITIAL  CONDITION,  ZO  xxxxxxxx 

B2=SQRT(UINFXDENSL/(DEFFXVISD) 

THETAC  0)  =  1 . E-9 

X(0)=1 .E-9 

CALL  DRV( F2, 0 . , F21 ) 

CALL  DRV( F3, 0 . , F31 ) 

Cl=(2XFXGl(0.)+DEFF/AXPXF21+3x(DEFF/A)xx2xFxB0XF31)/3. 

C2=F11 ( 0 . ) 

C3=-2 . XAXXFNC  0 . )/DEFF 

Z(0)=SQRT((SQRT(C2XX2-4XC1XC3)-C2)/(2.XC1)) 

ALPHAO=CONDLXB2/Z( 0 ) 

SUM=0. 

xxx  DETERMINE  DIMENSIONLESS  FILM  THICKNESS,  Z  xxxx 
xxx  ADAMS-MOULTON  PREDICTOR/CORRECTOR  ALGORITHM  xxxx 

FCO  =  FCN( THETAC  0 ) , Z( 0 ) ) 

ZPI=Z(0)+HXFC0 

Z(1)=Z(0)+H/2.X(FC0+FCN<THETA(1),ZPD) 

FC1  =  FCN( THETA( 1 ) , Z  C 1 ) ) 

ZP2=Z( 1 )+H/2 . x<  3 . XFC1-FCO ) 

Z(2)=Z(I)+H/I2.X( 5XFCNC THETAC 2) ,ZP2)+8 . XFC1-FC0 ) 

FC2=FCN( THETAC 2) ,Z(2) ) 

ZP3=ZC2)+H/I2.XC23.XFC2-I6.XFC1+5.XFC0) 

ZC  3) =ZC2)+H/24 . XC  9 . XFCNC THETAC 3) ,ZP3)+19 . XFC2-5 . XFCI+FCO) 
FC3=FCNC THETAC 3) , ZC3) ) 

DO  50, 1=4, M 

ZP4=ZC I-I >+H/24 . XC55 . XFC3-59 . XFC2+37 . XFC1-9 . XFCO ) 

ZC I )=ZC 1-1 )+H/720. XC 251 .XFCNC THETA Cl), ZP4)+646.XFC3-264.X 
+  FC2+I06 .XFC1-19 .XFCO) 

FC0=FC1 
FC1 =FC2 
FC2=FC3 

FC3=FCNC THETAC I ) ,Z(I ) ) 

C4=( 2 . XFXGC THETAC I ) )+DEFF/AxPXF2C THETAC I ) )+3 . X( DEFF/A)XX2XF 
+  XB0XF3C  THETAC I ) ) )xZC I ) XX3+F1 (THETAC I ) )XZ( I ) 

C5=C4+FI( THETAC I ) )XZ(I) 

IFCC4.LE.0.)  THEN 

PRINTX, 'C4  LIMIT' 

IC=I 

GO  TO  60 

ENDIF 

IFCC5.LE.0.)  THEN 

PRINTX, *C5  LIMIT* 

IC=I 

GO  TO  60 

ENDIF 
50  CONTINUE 
IC=N-1 

60  PRINTX, 'COMPLETED  Z  DETERMINATION' 

xxx  DETERMINE  HEAT  TRANSFER  PROPERTIES  xxxxxxxxx 


J  =  0 
JJ  =  0 

DO  100, 1=1, IC 
DELNDSZ( I ) 
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FILE:  MIXED1  FORTRAN  A 


DEL=DELND/B2 

ALPHA=CONDL/DEL 

SUM=SUM+! ALPHAOXAXXFN! THETA! 1-1 ) )+ALPHAXAXXFN!THETA! I ) ) )XHX . 5 

ALPHAO=ALPHA 

IF(JJ.EQ.J)  THEN 

QFLUX=ALPHAX(TS-TW) 

NU=DEFF/DEl 
DQFL  UX= 1 ,/DELND 

WRITE! 20, 5030)  THETA (I > *180 ./PI , X(I )X2 ./! DEFFXPI ) , 

+  DEL*1000 . , QFLUX/1 JOO . , DELND, DQFLUX, NU 

WRITE! 30,  JO  X(I)X2./(DEFFxPI), DELND 
J=J+1./1S0. XN+1 

ENDIF 
JJ=I+1 
100  CONTINUE 

WRITE! 20, 5030)  THETA!  10*180  ./PI , X!  10*2 ./(DEFFXPI ),  DEL*1 000  . , 
-t-QFLUX/1000  .  ,  DELND, DQFLUX, DEFF/DEL 
WRITE! 30,  SO  X! IC ) *2  ./! DEFFXPI ),  DELND 

NUAVG=SUMX2/(PIXC0NDL ) 

WRITE!20,5040)  ' NU/SQRT! RE)  =  ' . NUAVG/SQRT!RE) , 'NUAVG  =  ', NUAVG 

5000  FORMAT! IX, T5, 70! *X’ '  //  IX, T5, 'HEAT  TRANSFER  PROPERTIES'  //  2!1X, 
+T1 0 , A17 , F6 . 2  /  ),4!1X,T10,A17,E9.3  /)) 

5010  FORMAT!  /  IX, T5 , ' ELL  IPSE  GEOMETRY*  //  3(1X,T10,A19,F6 .4  /  )) 

£020  FORMAT!  /  IX, T7 ,' THETA ', T15 ,' DIM  X' , T25, ' DEL ' , T32, ' HT  FLUX’, 

+T41 , ' DIM  DEL',T52, T  DIM' , T58 , 'NUSS  *'  /  T7 , • ! DEG) • , T25, ’ !MM) ’ , T32, 
+ 1 ! KW/M2 ) ' , T50 , ’HT  FLUX'  /  T6 , 7 !'='), T15 , 6 !’='), T23, 7 !’=•), T32, 
+7!'  =  '),T4i  .7! '  =  ' ),T50,7! *  =  ' ),T58,6!  '  =  ')  /) 

5030  FORMAT !  1)  6  ,  F7 . 3 ,  T15 ,  F6 . 4, T23,  F7 . 4,  T32,  E8 . 2,  T41 ,  F7 . 4, T50 ,  F7 . 4 , 

+T58 , E9 . 3  ) 

5040  FORMAT!  //  3<1X, T5 , A17 , F8 . 4  /)) 

6000  FORMAT !1X,F4.1,F10.3,F7 .3) 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

FUNCTION  XFN( PHI ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  THIS  FUNCTION  SUBPROGRAM  IS  USED  TO  DETERMINE  STREAMWISE  LENGTH 

x  AS  A  FUNCTION  OF  THETA.  XFN  REPRESENTS  PARAMETRIC  RADIUS, 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  XFN,PHI,A,K,DEFF 
COMMON  /GEOM/A , K, DEFF 

XFN=SQRT! !KXX2+ . 25X! (KXX2-1 . )XSIN! 2XPHI ) )XX2)/( ! COS! PHI ) )XX2 
++! KXSIN! PHI ) )xx2) ) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  G(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxrxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  GRAVITY  FUNCTION, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  G, PHI, A, K, DEFF 
COMMON  /GEOM/A, K, DEFF 

G=SIN(PHI)/SQRT( (SIN! PHI) )XX2+( KXCOS! PHI ) )XX2) 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX^'XXXXX 
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FILE:  MIXED1 


FORTRAN  A 


FUNCTION  Gl(PHI) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

x  THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  STREAMWISE  DERIVATIVE 
x  OF  THE  GRAVITY  FUNCTION, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL  G1,PHI,A,K,D£FF 
COMMON/GEOM/A, K, DEFF 

G1=KXX2XC0S( PHI )/SQRT(((SIN(PHI) )XX2+(KXC0S(PHI ) )XX2)XX3) 

END 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXXXXXXX 

FUNCTION  FI (PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  VAPOR  VELOCITY 

x  FUNCTION 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  FI, PHI, A, K, DEFF 
COMMON  /GEOM/A,  K, DEFF 

F1=(1.+K)XSIN(PHI)/SQRT((SIN(PHI))XX2+(KXC0S(PHI))XX2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  Fll(PHI) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

x  THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  STREAMWISE  DERIVATIVE 
x  OF  THE  VAPOR  VELOCITY  FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  Fll , PHI. A, K, DEFF 
COMMON/GEOM/A, K, DEFF 

Fll=( 1 .+K)XKXX2XC0S (PHI) /SORT (((SIN(PHI) )XX2+(K*C0S(PHI) )XX2)*»3) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  F2( PHI ) 

XKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

x  THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  PRESSURE  GRADIENT 

X  FUNCTION. 

X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F2, PHI , A,K, DEFF 
EXTERNAL  XFN 
COMMON/ GEOM/A ,K, DEFF 

F2=((l .+K)XK)XX2XSIN(2XPHI)/((SIN(PHI))XX2+(KXC0S(PHI))XX2) 

+/XFN( PHI ) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  F3(PHI ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x  THIS  FUNCTION  SUBPROGRAM  REPRESENTS  THE  SURFACE  FUNCTION. 

X 

XXXXXXXXXXXXXKXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
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FILE:  MIXED1  FORTRAN  A 


REAL  F3,PHI,A,K,DEFF 
EXTERNAL  XFN 
COMMON/GEOM/A,K,DEFF 

F3=KX(1 .-KXX2)XSIN(2XPHIVSQRTCCCSIN(PHI))XX2+(KXC0SCPHI))XX2) 
+xx5)/XFN(PHI) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  FCN(PHI.Z) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  THIS  FUNCTION  SUBPROGRAM  IS  USED  IN  A  PROBLEM-SOLVER 
X  SUBROUTINE  TO  SOLVE  THE  O.D.E. 
x  DZ/DX  =  F(PHI.Z) 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  F21,F31,PHI,FCN1,FCN2,FCN,A,K,DEFF,F,B0,P 

COMMON  /GEOM/A, K, DEFF/PARAM/F, BO, P 

EXTERNAL  XFN,G,G1, FI , Fll, F2, F3 

CALL  DRV( F2, PHI , F21 ) 

CALL  DRV(F3,PHI,F31) 

FCN1=(2.XFXG(PHI)+DEFF/AXPXF2(PHI)+3.X(DEFF/A)XX2XFXB0XF3(PHI)) 

+XZXX3+F1(PHI)XZ 

FCN2=2 . XAXXFNC  PHI )/DEFF-( 2 . XFXG1 (PHI )+DEFF/ AXPXF21+3 . X( DEFF/A) xx2 
+XFXB0XF31 ) XZXX4/3 . -FI 1 ( PHI ) XZXX2 
FCN=FCN2/FCN1 

END 

XXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE  DRV(FN,PHI,FN1) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  THIS  SUBROUTINE  NUMERICALLY  DIFFERENTIATES  A  FUNCTION, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
REAL  FN, PHI , FN1 , EPS 

FN1=( FN(PHI+1 . E-6 ) -FN(PHI-1 . E-6 ) )/2 . E-6 

RETURN 

END 
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APPENDIX  C  -  COMPUTER  CODE  FOR  TWO-PHASE  BOUNDARY  LAYER 


MODEL 

*  PROGRAM  VAPBL1 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

x  THIS  PROGRAM  NUMERICALLY  SOLVES  THE  PROBLEM  OF  MIXED 
x  CONVECTION  CONDENSATION  WITH  SURFACE  TENSION  ON  AN  ELLIPTICAL 
X  CYLINDER.  THE  PROBLEM  FORMULATION  USES  THE  SOLUTION  OF  THE 
X  TWO-PHASE  BOUNDARY  LAYER  EQUATIONS  FOR  THE  CONDENSATE  AND  VAPOR. 
X 

x  THE  PROGRAM  USES  AN  IMSL  NUMERICAL  PROBLEM  SOLVER,  DIVPAG, 
x  TO  SOLVE  A  1ST  ORDER  ORDINARY  DIFFERENTIAL  EQUATION.  THE 
x  PROBLEM  SOLVER  IS  WRITTEN  IN  DOUBLE  PRECISION  WHILE  THE  MAIN 
x  PROGRAM  IS  IN  SINGLE  PRECISION.  THE  MAIN  PROGRAM  MUST  BE 
x  COMPILED  IN  DOUBLE  PRECISION  (AUTO  DOUBLE) . 
x 

x  THIS  PROGRAM  CAN  BE  RUN  WITH  DIMENSIONLESS  PARAMETERS 
x  ( F,G)  SPECIFIED  OR  WHERE  THE  PARAMETERS  ARE  CALCULATED 

x  BASED  ON  SPECIFIED  FLUID  PROPERTIES.  MODIFICATIONS  TO  THE 

x  PROGRAM  (REMOVAL  OF  UNNECCESSARY  STEPS  BY  "COMMENTING"  THE 

X  STATEMENT  OUT  OF  THE  PROGRAM)  IS  REQUIRED.  SPECIFIED  FLUID 

x  PROPERTIES  INCLUDE  T  SAT,  T  WALL,  AND  UINF. 
x 

X  IF  THE  ALGORITHM  HAS  DIFFICULTY  LOCATING  THE  VAPOR  SEPARATION 

x  POINT,  THEN  THE  PROGRAM  NEEDS  TO  BE  MODIFIED  TO  INPUT  THE 

X  PARAMETRIC  AND,  PHI,  FOR  THE  SEPARATION  POINT, 
x 

x  NUMERICAL  RESULTS  ARE  DIRECTED  TO  AN  EXTERNAL  DATA  FILE. 


x 

x  MAJOR  VARIABLES  USED  ARE: 
x 


X 

PHI 

PARAMETRIC  ANGLE 

X 

X 

STREAMWISE  LENGTH 

X 

DEFF 

EFFECTIVE  DIAMETER 

X 

DEL 

CONDENSATE  FILM  THICKNESS 

X 

DELND 

DIMENSIONLESS  FILM  THICKNESS 

X 

z 

DIMENSIONLESS  VAPOR  MOMENTUM  THICKNESS 

X 

F, JARPRL 

DIMENSIONLESS  PARAMETERS 

X 

UINF 

VAPOR  FREE  STREAM  VELOCITY 

X 

TS 

VAPOR  SATURATION  TEMPERATURE 

X 

TW 

WALL  SURFACE  TEMPERATURE 

X 

DENSL 

CONDENSATE  DENSITY 

X 

DENSV 

VAPOR  DENSITY 

X 

CONDL 

CONDENSATE  THERMAL  CONDUCTIVITY 

X 

VISL 

CONDENSATE  VISCOSITY 

X 

VISV 

VISCOSITY  OF  VAPOR 

X 

HFG 

LATENT  HEAT  OF  VAPORIZATION 

X 

NU 

NUSSELT  NUMBER 

X 

REV 

VAPOR  REYNOLDS  NUMBER 

X 

RE 

TWO-PHASE  REYNOLDS  NUMBER 

X 

NUAVG 

MEAN  NU  NUMBER 

X 

ALPHA 

HEAT  TRANSFER  COEFFICIENT 

X 

QFLUX 

HEAT  FLUX 

X 

A 

SEMI-MAJOR  AXIS 

X 

K 

ECCENTRICITY 

X 

XFN 

PARAMETRIC  RADIUS 

X 

G 

GRAVITY  FUNCTION 

X 

FI 

VELOCITY  FUNCTION 

X 

Fll 

DERIVATIVE,  VELOCITY  FUNCTION 

X 

F2 

PRESSURE  GRADIENT  FUNCTION 

X 

F3 

SURFACE  TENSION  FUNCTION 

X 

TAU 

INTERFACIAL  SHEAR 

X 

KAPPA 

PRESSURE  GRADIENT  FORM  PARAMETER 

X 

KAPPAA 

SUCTION  FORM  PARAMETER 

X 

XXXXXXXXXKXXXXXXXXXXXXXXKXKKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

PARAMETER! PI *3 . 141592654, GRAV»9 . 807 , NEQ*2, NEQ1=1 , NPARAM*50 ) 

REAL  PHI(0:180000),DEL(0: 180000 ),x(0: 180000 ),A,TS,TW, 
+H,A1,AZ,A3,A4,A5,A6.A7,A8,A9,A10,A11,TF, PS AT ,TB,TC,TD,R,TE, DENSV , 
+TL , DENSL , VISL , VISV , TG, CONDL , HFG, SIGML ,F,K,X1,X2,X3,DEFF,UINF, 

♦REV , RE, JARPRL , DELOO , DEL  0 ,Z0 , TAU,  DTAU, FDELO , DFDELO , DELI ,  Z1 , ALPHAO, 
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nn  on 


FILE:  VAPBL 


FORTRAN  A 


+SUM,AA(1,1) , FCN,FCNJ, HINIT,PARAMC NPARAM), KAPPA, KAPPAA, DKAPPA, 

+T0L , XEND, XBEG, Y(2) , DDEL , DZ, Yl( 1 ) , NU, NUAVG, NUORE, DKAPPAA , XO, DXO , 
+DELND, DELL , ALPHA, QFLUX, DQFLUX, D2UPHI ,  DUPHI ,  UPHI 

INTEGER  N»  JJ , II, IDO , IMETH 

COMMON  /GEOM/A , K, DEFF/PARAM/ JARPRL , F 

EXTERNAL  ZFN, DELFN, XFN, G, G1 , FI , Fll , NUSSELT, DRV, FCN, FCNJ , IVPAG, 
+SSET 

QPEN(20, FILE3 '/VAPOR1  OUTPUT  A1 ', STATUS3 ' OLD' ) 

OPEN ( 30 , FILE3 '/VAP0R2  OUTPUT  A1 ', STATUS3 'OLD' ) 

PRINT*, ' ELLIPSE  OF  THE  FORM  (X/B )**2+(Y/A)**23l » 

PRINT*, 'INPUT  "A"  DIMENSION  (M)T' 

READ*, A 

PRINT*, 'INPUT  T  SAT  AND  T  WALL  CIO?’ 

READ*,TS,TW 

PRINT*, 'INPUT  ANGLE  STEP  SIZE  (DEGREES)?' 

READ*,H 

PRINT*, 'STEP  SIZE,  H  3  ',H 
PRINT*, ’INPUT  PHI  AT  SEPARATION?' 

READ*, PHIS 

PRINT*, ’INPUT  UINF?' 

READ*, UINF 
PRINT*, 'INPUT  <?• 

READ*, K 

PRINT*, 'INPUT  G?' 

READ*, JARPRL 

N=180/H+l 
PRINT*, 'N-1=’,N-1 
H=PI/180.*H 

***  DETERMINE  FLUID  PROPERTIES  *x*xxxx* 

Al=15. 49217901 
A2=-5. 6783717693 
A3=l. 4597584637 
A4=13. 877000608 
A5=-80. 887673591 
A6=123. 56883468 
A7=-188. 321212064 
A8=660. 91763485 
A9=-1382. 4740091 
A10s1300 .1040184 
All=-449. 39571976 
TF3TS/1000. 

PSAT=1 . E+6*EXP( A1+A2/TF+A3*L0G(TF)+TF*A4+A5*TF*X2+A6XTF**3+A7X 
+TF**4+A8XTF**5+A9*TF*X6+A10XTF**7+A11*TF*X8) 

TB31500 ,/TS 

TC32 . 5*L0G( 1 . -EXP( -TB) ) 

TD3.0015/(1.+.0001xTS)~. 000942*EXP(TB+TC)*SQRT( 1 ./TB)- . 000488 2*TB 
R3461 .51 

TE=2 . *PSAT/(RXTS ) 

DENSV=TE/C1 . +SQRTC1 .+2.XTDXTE) ) 

TL=(2.*TW+TS)/3. 

DENSL  =  1 ./( . 0012674-TLXC  2 . 02915E-6-TLX3 .8333E-9) ) 

VISL32 . 414E~5*10XX( 247 . 8/CTL-140.  )) 

VISV=-4. 478415E-6+TSXC5 . 0216E-8-1 .579E-11XTS) 

TG=TL/273 . 15 

CONDL3- . 92407+TG*(2.8395-TG*( 1 .8007-TGx( . 52577-TGX. 07344) ) ) 
HFG=3468920 . -TSX( 5707 . 4-TS*( 11 . 5562-TSx .0133103)) 

SIGML=(- . 0003*(TL-273 . 15)X*2- . 138XCTL-273 . 15)+75 .6 )/1000 . 

**x  DETERMINE  PHK I ) , X( I ) , DEFF  xxxxxxxx 

*xx  USES  SIMPSON'S  RULE  TO  NUMERICALLY  INTEGRATE  ALONG  SURFACE  *xx 

PHI(0)=0. 

X(0)30. 
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FILE:  VAP8L 


FORTRAN  A 


DO  10,  1=1, N 

PHI(I)=PHI(I-1)+H 

X1=XFN(PHI(I-1)) 

X2=XFN(  (PHICI-l  )+PHI(I) )X.5) 

X3=XFN(PHI< I ) ) 

X( I ) =X( 1-1 )+H/6 . XCX1+4.XX2+X3) 

10  CONTINUE 

DEFF=2.XXCN)/PI 

PRINTX, 'COMPLETED  DE  DETERMINATION,  DE=  * , DEFF 

xxxx  DIMENSIONLESS  PARAMETER  DETERMINATION  xxxx 

F=VISLXDEFFXHFG/(C0NDLXUINFXX2*(TS-TW)) 

REV=UINFXDENSVXDEFF/VISV 

RE=UINFXDENSLXDEFF/VISL 

JARPRL=CONDLX( TS-TW)/HFGXSQRT ( DENSL/C  DENSVXVISLXVISV) ) 

PRINTX, 'COMPLETED  PARAMETER  DETERMINATION' 

WRITE(20,5000)  'SAT.  TEMP.  CIO  =  ',TS,'WALL  TEMP.  =  ',TW, 

+'F  =  ' , F, ' UINF  (M/S)  =  ' ,UINF, 'G  =  *,JARPRL 
WRITEC20, 5010)  '"A"  LENGTH  CM)  =  ',A,",IC  ECCENTRICITY  =  ',K, 
+'DEFF  CM)  =  ' , DEFF 

xxx  DETERMINE  INITIAL  VALUES  OF  Z  AND  DEL  xxx 

XXX  USES  INTIAL  VALUE  OF  DEL  FROM  SHEKRILADZE-GOMELAURI  xxx 

xxx  FOR  MIXED  CONVECTION  AS  STARTING  POINT  TO  FIND  Z  THEN  xxx 

xxx  REDETERMINES  DEL  FOR  THE  PREDICTED  VALUE  OF  Z  USING  xxx 

xxx  NEWTON'S  METHOD,  CHECKS  FOR  CONVERGENCE  AND  ITERATES  xxx 

DEL00=SQRT( . 75XK/FX(S<5RT( (Cl. +K)/K)xx2+16 . xFx A/ ( 3 . XDEFF) ) 
+-C1.+0/K)) 

DEL0=DEL00 

Z0=C-.0735XXFN<0. )XJARPRL/DEFF/F11C0.)/DEL0+SQRT(( .0735X 
+XFNC0. )XJARPRL/DEFF/F11C0.)/DEL0)XX2+.0735XXFNC0.)/DEFF/ 

+F11C0. )))*X2 

UPHI=F1(0 . ) 

DUPHI=F11 ( 0 . ) 

CALL  DRVC  Fll , 0 . , D2UPHI) 

X0=XFN(  0 . ) 

CALL  DRVCXFN, 0 . , DXO) 

20  KAPPA=DEFF/X0XDUPHIXZ0 

30  KAPPAA=.0«82+.17AXJARPRLXSQRTCZ0)/DEL0 

DKAPPA=DEFFXZOXD2UPHI/XO-DEFFXDUPHIXZOXDXO/XOXX2 

DTAU=6.44XSGRTCKAPPAXKAPPAA+KAPPAAXX2)XDUPHI/S0RTCZ0)+3.22XUPHI 

+XKAPPAAXDKAPPA/SQRT(Z0)/SQRTCKAPPAXKAPPAA+KAPPAAXX2) 

FDEL0=X0/DEFF-F/3.XG1(0. )XDEL0XX4-DTAUXDEL0XX3/C4.XJARPRL) 
DFDEL0=-4.XFXG1(0 . )XDEL0XX3/3.-3.xDTAUxDElQxx2/C4.xjARPRL) 
DEL1=DEL0-FDEL0/DFDEL0 

IFC CABS (DELI -DEL 0) .GT. < . OOlxDELQ ) ) .AND. CFDELO . GT. 1 . E-7 ) )  THEN 
DEL0=DEL1 
GO  TO  30 
ENDIF 

Z1=<-.0735XXFN<0. )XJARPRL/CDEFFXF11C0.)XDEL1)+SQRTCC .0735X 
+XFN(0.)XJARPRL/(DEFFXF11(0. )XDEL1))XX2+.0735XXFN(0. )/(DEFFX 
+F11 C  0 . ) ) ) )XX2 

IF  C CABS C DELI -DEL 00) ,GT. ( . OOlxDELOO)) .OR. (ABSCZ1-Z0) . GT. C . OOlx 
+Z0)))  THEN 
Z0  =  Z1 

DEL00=DEL1 
DEL0=DEL1 
GO  TO  20 
ENDIF 

DELC0)»DEL1 
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FILE:  VAPBL 


FORTRAN  A 


ALPHAO  =  l  ./DEL(O) 

SUM=0. 

PRINT*, 'COMPLETED  ZCO),DELIO)  DETERMINATION*,'  Z=',Z1, 
+'DEL=',DEL<0) 

XXX  EVALUATION  OF  Z(I)  AND  DEL(I)  ***x 

xxx  USES  IMSL  ROUTINE  DIVPAG  xxxx 

HINIT=l.£-4 

IMETH=2 

CALL  SSETCNPARAM, 0.0, PARAM, 1 ) 

PARAMC 1 ) =HINIT 

PARAM(12)=IMETH 

XBEG=PHI(0) 

YC1)=Z1 

Y(2)=D£L(0) 

TOL=l . E-3 
IDO=l 

DO  100,1=1, N-l 

C  IFCPHI(I) .GT.PHIS)  GO  TO  190 

XEND=PHI C 11 

CALL  DIVPAGL IDO , NEQ, FCN, FCNJ , AA, XBEG, XEND, TOL , PARAM, Y) 
DEL(I)=Y(2) 

KAPPA=D£FF/XFN<PHI(I))*F11(PHI(I))*YC1) 

KAPPAA= . 0682+ . 174*JARPRL*SQRTC Y( 1 ) )/YC2) 

TAU=6 . 44XF1 C  PHI C 1 3  3 XSQRTCKAPPAXKAPPAA+KAPPAAX*2)/SQRT( Y( 1 ) ) 

WRITEC30,*)  'PHI=',PHICI)X180/PI, 'TAU=',TAU 

II=I 

xx*  CHECK  FOR  VAPOR  BOUNDARY  LAYER  SEPARATION  xxxx 

IFCTAU.LE.O.)  GO  TO  190 
100  CONTINUE 

190  PRINTX, 'SEPARATION  ANGLE  AT  PHI  = ’ , PHIC II)*180 ./PI 
PRINT*, «TAU= ' , TAU, ' DEL  = ' , DEL C II) 

ID0=3 

CALL  DIVPAGCIDO, NEQ,FCN, FCNJ, AA, XBEG, XEND, TOL, PARAM, Y) 

XXX  USES  A  SIMPLE  NUSSELT  ANALYSIS  TO  COMPUTE  FILM  THICKNESS  **** 
XXX  DOWNSTREAM  OF  THE  SEPARATION  POINT.  PROBLEM  SOLVER  IS  XXXX 

xxx  IMSL  DIVPAG  xxxx 

IFCII . HE. (N-l) )  THEN 
XBEG=PHI(II) 

IDO=l 

Yl< 1)SDEL(II) 

DO  200,  1=11+1, N-l 
XEND=PHI ( I ) 

CALL  DIVPAGL IDO,  NEQl, NUSSELT, FCNJ, AA, XBEG, XEND, TOL, PARAM 
+  ,Y1) 

DEL(I)=Y1<1) 

200  CONTINUE 

ENDIF 
ID0=3 

CALL  DIVPAGC IDO , NEQ1 , NUSSELT , FCNJ , AA, XBEG, XEND, TOL , PARAM, Y1 ) 

PRINTX, 'COMPLETED  ZCI),  DELCI)  DETERMINATION' 

WRITEC20, 5050)  PHI<II)*180 ./PI 
WRITEC20, 5020) 

xxx  DETERMINE  HEAT  TRANSFER  PROPERTIES  xxxxxxxxx 


JJ=0 

DO  300, 1=1, N-l 
DELND=DELCI) 

DELL=DELNDXS0RT(VISLXDEFF/(DENSLXUINF)) 

ALPHA=1./DEL(I) 

SUM=SUM+(ALPHA0XXFN(PHI(I-1))+ALPHAxXFN(PHI(I)))xHX.5 

ALPHAO=ALPHA 


FILE:  VAPBL 


FORTRAN  A 


IF(CJJ.EQ.I).OR.CJJ.EQ.O))  THEN 
QFLUX=ALPHAX(TS-TW) 

NU=DEFF/DELL 
DQFLUX=1 ./DELND 
JJ=JJ+10./130.XN+1 

WRITEC2Q , 5030 )  PHI ( I )X180/PI , XC I )X2/(DEFFxPI) , DELL 
+  XI 000, QFLUX/1000, DELND, DQFLUX,NU 

ENDIF 

300  CONTINUE 

WRIT£(20,5Q30)  PHI(N-l)xl80./PI, XCN-l )*2./( DEFFXPI ) , DELLX1000 . , 
+QFLUX/1Q00 . , DELND, DOFLUX, DEFF/ DELL 
WRITEC30,*)  X(N-l)X2./( DEFFXPI ) , DELND 

NUAVG=SUMX2/(PIXDEFF) 

WRITE(20,5040)  ’NU/SQRT(RE)  =  ' , NUAVG/SQRT ( RE) 

5000  FORMAT ( IX, T5,70( '»•)  //  IX, T5, 'HEAT  TRANSFER  PROPERTIES’  //  2(1X, 
+T10,A17,F6 .2  /  )»5(1X,T10,A17,E9.3  /)) 

5010  FORMAT (  /  IX, T5 , ’ ELL  IPSE  GEOMETRY’  //  3C1X.T10, A19 , F6 . 4  /  )) 

5020  FORMAT C  /  IX, T7 , ’THETA ’, T15, ’ DIM  X • , T25, • DEL ’ , T32, 'HT  FLUX’, 

+T41 , ’ DIM  DEL ’,T52, 1  DIM* , T58, ' NUSS  *'  /  T7 , ’ ( DEG) ’ , T25, ’ (MM) ’ , T32, 
+’<KW/M2)’,T50,'HT  FLUX’  /  T6 , 7 (’='), T15, 6( ’=') ,T23, 7 (’='), T32, 

+7C ’=’),T41,7(’=' ),T50,7C*a* ),T58,6< ’=*)  /) 

5030  FORMAT! IX, T6 , F7 . 3,T15, F6 .4, T23, F7.4,T32,E8.2,T41,F7.4,T50,F7.4, 
+T58, F9 . 1 ) 

5040  FORMA T(  //  3C1X, T5, A17 , F8 . 4  /)) 

5050  F0RMAT(1X,T6, 'VAPOR  SEPARATION  AT  PHI  =  ’,F6.2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  XFN(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  THIS  FUNCTION  SUBPROGRAM  IS  USED  TO  DETERMINE  STREAMWISE  LENGTH 
X  AS  A  FUNCTION  OF  PHI. 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  XFN, PHI , A, K, DEFF 
COMMON  /GEOM/A, K, DEFF 

XFN=AXSQRT((KXX2+.25X((KXX2-1.)XSINC2XPHI))XX2)/((C0SCPHI))XX2 
++(KxSIN(PHI) )**Z) ) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  GCPHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x  GRAVITY  FUNCTION 
X 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAL  G, PHI, A, K, DEFF 
COMMON  /GEOM/A, K, DEFF 

G=SIN(PHI )/SQRT( ( SIN(PHI) )XX2+(K*C0S(PHI ) )XX2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  Gl(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  STREAMWISE  DERIVATIVE  OF  GRAVITY  FUNCTION 
x 
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FILE:  VAPBL 


FORTRAN  A 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

REAL  G1 , PHI , A, K, DEFF 
COMMON/GEOM/A, K, DEFF 

G1=KXX2XC0S(PHI )/ SORT ((<SIN( PHI) )XX2+(KXC0S(PHI ) )XX2)XX3) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  Fl(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  VAPOR  VELOCITY  FUNCTION 
X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  FI, PHI, A,K, DEFF 
COMMON  /GEOM/A, K, DEFF 

Fl=<l.+K)xSIN<PHI)/SQRT<CSIN(PHI))xx2+CKxC0S(PHI))xx2) 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  Fll(PHI) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x  STREAMWISE  DERIVATIVE  OF  VELOCITY  FUNCTION 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  Fll , PHI , A, K, DEFF 
COMMON/GEOM/A , K , DEFF 

Fll  =  a.+K)XKXX2XC0S(PHI)/S9RT(((SIN(PHI))XX2+(ICXC0S(PHI))XX2)XX3) 
END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
FUNCTION  DELFN(PHI,Z,DEL,DZ) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  STREAMWISE  DERIVATIVE  OF  CONDENSATE  FILM  THICKNESS 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  DEL FN, PHI , Z, DEL , DZ, DELFNI , DELFNZ, DELFN3, X, FF1 , FF11 , FF12, DXFN 
+, A , K, DEFF, J ARPRL , F, KAPPA , KAPPAA , XX, FKAPPA 

EXTERNAL  XFN, FI , Fll ,G,G1 


COMMON/GEOM/A, K, DEFF/PARAM/JARPRL , F 

IFCPHI . EQ . 0 . D+O)  THEN 
DELFN=0 . D+Q 
ELSE 

FF1=F1CPHI ) 

FF11=F11(PHI) 

CALL  DRV ( Fll , PHI , FF12) 

X=XFN( PHI ) 

CALL  DRV(XFN, PHI, DXFN) 

KAPPA=DEFFXFF11XZ/X 

KAPPAA* . 0682+ . 174XJ ARPRLXSQRT(Z)/DEL 

XX=KAPPAAXKAPPA+KAPPAAXX2 

FKAPPA=$QRTCXX) 

DELFNI*. 805XKAPPAAXDEFFXFFlxFFllXDELXX3/(JARPRLxFKAPPAXXXSqRT(Z)) 

++(7 .0035E-2)X(2.XKAPPAA+KAPPA)XFF1XDELXX2/(FKAPPAXZ)-.805XFKAPPAX 

+FF1xDEL**3/( JARPRLXZXX1 .5) 
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FILE:  VAPBL 


FORTRAN  A 


DELFN2=FXGCPHI  )#DEL**3+3 . 22XFFlXFKAPPAxDELXX2/( JARPRLXSQRT(Z) ) 

+-.  14007X(2.XKAPPAA+KAPPA)XFF1XDEL/FKAPPA 

D£LFN3=X/DEFF-FXG1 (PHI)* DEL **4/ 3.-1. 6 1XFKAPPAxFF11xDELXx3/( JARP ; L 
+XSQRTCZ) )- .305XKAPPAAXDEFFXFF1XFF12XSQRT<Z)XDELXX3/C JARPRLXFKAPPAX 
+X)+.305XKAPPAAxDEFFxDXFNxFFlxFFllxsqRT(Z)XDELxx3/C JARPRLXFKAPPAX 
+XXX2) 

DELFN=(DELFN3-DELFN1XDZ)/DELFN2 

ENDIF 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

FUNCTION  ZFN<PHI,Z,DEL) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX'  ‘XXXXXXXXXXXXXXXXXXXXXXXXXXX 

X 

X  STREAMWISE  DERIVATIVE  OF  VAPOR  MOiW.rUM  THICKNESS  PARAMETER 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  ZFN, PHI ,Z,  DEL , DEFF, A, K, JARPRL ,  F 

EXTERNAL  XFN 
EXTERNAL  FI 
EXTERNAL  Fll 

COMMON/GEOM/A, K, DEFF/PARAM/JARPRL, F 

IFCPHI . EQ . 0 . D+0)  THEN 
ZFN=0 . D+O 
ELSE 

ZFM=C .441x<1.-2.XJARPRLXSQRT(Z)/DEU-6.XDEFFxF11(PHI)xZ/XFN(PHI)) 
+XXFNCPHI)/F1(PHIVDEFF 
ENDIF 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE  NUSSELT(NEQ,PHI,Y,YPRIME) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  STREAMWISE  DERIVATIVE  OF  CONDENSATE  FILM  THICKNESS  USING 
X  A  NUSSELT  ANALYSIS  DOWNSTREAM  OF  VAPOR  SEPARATION  POINT, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

INTEGER  NEQ 

REAL  PHI , Y(NEQ) , YPRIMEC NEQ) , A, K, DEFF, JARPRL , F 
EXTERNAL  XFN,G,GI 

COMMON/GEOM/A, K, DEFF/PARAM/JARPRL, F 

YPRIMEC1)=(XFNCPHI)/DEFF-F/3.XYC1)*X4XG1CPHI))/(FXG(PHI)XYC1)XX3) 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE  DRV<FN,PHI,FN1) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

x  SUBROUTINE  TO  NUMERICALLY  DIFFERENTIATE  A  FUNCTION 

x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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FILE:  VAPBL 


FORTRAN  A 


REAL  FN, PHI , FN1 , EPS 

FN1 =(FN(PHI+l.E-6)-FN< PHI-1 ,E-6))/2.E-6 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE  TAUFNCPHI ,Z, DEL, DZ, DDEL , TAU, DTAU) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  SUBROUTINE  TO  DETERMINE  THE  INTERFACIAL  SHEAR  STRESS  AND 
x  ITS  STREAMWISE  DERIVATIVE, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

REAL  PHI , Z, DEL , DXFN , DZ, DDEL , KAPPA , KAPPA1 , KAPPAA, DKAPPA , 

+DKAPPA1 , DKAPPAA, TAU, DTAU, A, K, DEFF, JARPRL , F 

EXTERNAL  XFN,F1,F11 

CUMMON/GEOM/A, K, DEFF/PARAM/ JARPRL, F 

CALL  DRV( Fll , PHI , FI 21 
CALL  DRV(XFN, PHI, DXFN) 

KAPPA=DEFFXFI1(PHI)XZ/XFN(PHI) 

KAPPA1 = JARPRL XSQRT(Z)/ DEL 
KAPPAA= . 0682+ . 17  4XKAPPA1 

DKAPPA=DEFF/XFNCPHI)x(ZxF12+F11CPHI)xDZ)-DEFF/(XFN(PHI))XX2X 
+F11CPHI )XZXDXFN 

DKAPPA1=JARPRLx(.5xDELXDZ-ZXDDEL)/(SQRT(Z)xDELXX2) 

DKAPFAA= . 174XDKAPPA1 

TAU=6.44XF11(PHI)XSQRT(<KAPPAAXKAPPA+KAPPAAXX2)/Z) 

DTAU=6 . 44X( ( ( KAPPAAXKAPPA+KAPPAAXX2)XF11(PHI)+ . 5XFHPHI )x 
+(KAPPAAXDKAPPA+(2.XKAPPAA+KAPPA)XDKAPPAA))/SQRTCZX(KAPPAAX 
+KAPPA+KAPPAAXX2))-.5X(KAPPAAXKAPPA+KAPPAAXX2)XF1(PHI)XDZ/ 
+SQRT(ZXX3X(KAPPAAXKAPPA+KAPPAAXX2))) 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE  FCN(NEQ,X, Y,YPRIME) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X 

X  SUBROUTINE  TO  INPUT  DIFFERENTIAL  EQUATION  TO  BE  SOLVED 
X  BY  THE  PROBLEM  SOLVER, 
x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTEGER  NEQ 

REAL  X,Y(NEQ) .YPRIMECNEQ) 

EXTERNAL  ZFN, DELFN 

YPRIME( 1 )=ZFN(X, Y( 1 ) , Y<2) ) 

YPRIME(2) SDELFN(X, Y( 1 ) , Y( 2) , YPRIMEC 1) ) 

RETURN 

END 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
SUBROUTINE  FCNJ(NEQ,X, Y,DYPDY) 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
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*  tt  tt 


FILE:  VAPBL 


FORTRAN  A 


SUBROUTINE  REQUIRED  BY  PROBLEM  SOLVER,  NO  FUNCTION  IN  ANALYSIS 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
INTEGER  NEQ 

REAL  X, Y(NEQ) , DYPDY(X) 

RETURN 

END 
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